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SUMMARY 


We  have  investigated  the  potential  of  photoreffactive  beam  fanning  as  an  optical 
limiting  device.  Beam  fanning  is  attractive  as  an  optical  limiter  since  it  provides  the  possibility 
of  allowing  the  outside  world  to  be  monitored,  while  it  simultaneously  rejects  a  laser  threat. 
As  a  device  it  also  has  the  potential  to  be  sensitive  over  a  broad  spectral  range,  operate  with  a 
large  field-of-view,  provide  protection  against  continuous  and  pulsed  laser  sources,  protect 
against  a  laser  threat  from  any  direction,  simultaneously  block  multi-wavelength  laser 
radiation,  and  to  transmit  a  scene  without  significant  degradation  or  attenuation. 

The  major  obstacle  preventing  the  beam  fanning  limiter  from  reaching  this  potential 
has  been  speed  of  response.  Because  the  photoreffactive  beam  fanning  phenomenon  involves 
charge  transport,  it  is  inherently  slow.  Its  relatively  slow  speed  and  the  corresponding 
reputation  it  has  acquired,  suggests  that  the  beam  fanning  limiter  may  not  switch  fast  enough 
to  stop  a  laser  threat  before  sensor  damage  occurs. 

In  the  present  project  we  have  developed  several  techniques  which  are  successful  in 
enhancing  the  photoreffactive  beam  fanning  response  time.  In  particular,  we  demonstrate  and 
characterize  the  use  of  applied  d.c.and  a.c.  electric  fields,  use  of  a  grating  to  seed  the  beam 
fan,  use  of  cylindrical  and  spherical  lens,  control  of  the  acceptor  and  donor  concentrations, 
and  control  of  crystal  temperature.  Each  can  lead  to  an  improved  response  time. 

As  part  of  the  project  we  carried  out  a  demonstration  of  the  beam  fanning  optical 
limiter.  Using  a  combination  of  spherical  and  cylindrical  lenses  and  an  applied  d.c.electric 
field,  we  demonstrated  the  beam  fanning  limiter  for  a  laser  threat  of  lm  w/cm2  over  a  two  inch 
diameter  aperture.  The  limiter  switched  in  6  msec  with  an  O.D.  of  3  to  4.  This  limiter, 
therefore,  allows  only  6  ujoules  to  be  transmitted  before  the  O.D.  of  3  to  4  was  in  effect.  For 
both  higher  and  lower  intensities  the  limiter  switched  faster  and  slower,  as  expected,  still 
limiting  at  6  ujoules  with  an  O.D.  of  3  to  4.  The  demonstration  was  at  the  488  nm  and  514.5 
nm  wavelengths.  While  the  limiter  operated,  the  outside  world  could  still  be  monitored,  a 
large  field-of-view  was  in  effect,  and  about  50%  absorption  existed  for  outside  world  scenes. 


TABLE  OF  CONTENTS 


I .  OVERVIEW  OF  THE  PROBLEM  1 

II.  PHOTOREFRACTION  8 

III.  INCREASING  THE  INTENSITY  18 

IV.  APPLIED  DC  ELECTRIC  FIELD  26 

V .  APPLIED  A.C.  ELECTRIC  FIELD  6 1 

VI.  DETERMINATION  OF  TRAP  NUMBER  DENSITY  89 

VII.  GRATING  WRITING  AND  ERASURE  TIMES  109 

VIII.  SELF- FOCUSING  AND  DEFOCUSING  129 

IX.  CONTROL  OF  NA  AND  ND.  140 

X.  SEEDING  THE  BEAM  FAN  145 

XI.  SUMMARY  153 

XII.  PUBLICATIONS  (1990  -  TO  DATE)  158 

XIII.  PARTICIPATING  SCIENTIFIC  PERSONNEL  160 

XIV.  REPORT  OF  INVENTIONS  161 

XV.  BIBLIOGRAPHY  162 


I. 


OVERVIEW  OF  THE  PROBLEM 


The  existence  of  high  powered  lasers  is  a  severe  threat  to  all  optical  sensors. 
Widespread  use  of  lasers  as  range  finders  and  target  designators  makes  this  threat  a 
present  and  immediate  danger.  The  human  eye,  as  well  as  man-made  detectors  such  as 
the  many  FUR  devices  presently  employed  in  the  battlefield,  are  extremely  vulnerable 
to  stray  and  directed  laser  radiation.  The  eye  is  easily  attacked  by  light  from  a  laser 
source  at  a  considerable  distance.  Damage  caused  to  optical  sensors  is  extensive  and 
crippling.  With  the  success  of  many  military  operations  dependent  on  effective  use  of 
optical  sensors,  the  development  of  sensor  protection  is  imperative. 

The  laser  threat  itself  is  highly  agile.  The  high  powered  output  radiation  can 
occur  in  a  short  pulse  or  continuously,  while  the  wavelength  of  the  output  radiation  can 
be  either  fixed,  shifted,  or  mixed.  One  example  of  a  threat  laser  is  the  short-pulsed 
neodymium-doped  yittrium  aluminum  garnet  or  YAG  laser.  The  output  wavelength  of 
the  laser  is  in  the  near  infrared  at  about  one  micron.  However,  it  can  easily  be  altered 
via  frequency  doubling,  a  process  which  produces  an  additional  output  in  the  green  at 
about  0.5  microns.  The  output  can  be  further  modified,  producing  laser  radiation  at  a 
series  of  wavelengths  scattered  over  the  visible  spectrum  via  the  process  called  Raman 
scattering.  This  type  of  laser  threat  can,  therefore,  deliver  several  laser  frequencies 
simultaneously  onto  a  sensor.  It  is  commonly  referred  to  as  a  “rainbow”  or  “white- 
light”  laser.  The  actual  spectrum  of  the  output  laser  light  can  be  varied  and  difficult  to 
predict. 

While  the  output  of  the  YAG  laser  is  usually  a  short  pulse,  there  are  also  serious 
laser  threats  which  can  provide  a  continuous  output  of  laser  light.  The  argon-ion  and 
krypton-ion  lasers  are  prime  examples.  These  lasers  emit  radiation  at  seven  to  eight 
wavelengths  simultaneously  and  can  be  regarded  as  continuous  wave  “rainbow”  or 
“white-light”  lasers.  While  the  spectral  output  of  these  lasers  is  fixed,  the  development 


1 


of  continuously  tunable  dye  lasers  and  solid  state  lasers  presents  a  threat  with  a 
complex  and  variable  spectral  structure. 

The  ideal  protector  must  provide  complete  blocking  of  laser  radiation  above  a 
predetermined  threshold  value.  It  must: 

(i)  be  sensitive  over  a  broad  spectral  range; 

(ii)  not  significantly  degrade  or  attenuate  the  desired  scene  to  be 
interrogated; 

(Hi)  operate  with  a  large  field-of- view; 

(iv)  provide  protection  against  continuous  and  pulsed  laser 
sources; 

(v)  protect  against  laser  radiation  coming  from  any  direction; 

(vi)  be  capable  of  simultaneously  blocking  multi-wavelength  laser 
radiation  emitted  from  a  single  source. 

No  such  protection  device  presently  exists. 

Protection  devices  or  “optical  limiters”  can  be  classified  as  one  of  two  types: 
active  or  passive.  Active  optical  limiters  require  an  auxiliary  feedback  system  which  is 
initiated  by  the  incoming  laser  radiation.  A  simple  example  of  an  active  optical  limiter  is 
a  mechanical  shutter  which  closes  when  incident  laser  radiation  is  detected  by  an  optical 
diode  and  a  corresponding  electrical  signal  is  sent  to  the  shutter.  Since  lasers  with 
nanosecond  rise  times  are  now  common,  there  is  very  little  time  for  an  active  protection 
device  to  switch.  However,  electro-optic  switches  make  active  switches  a  realistic 
possibility.  Meanwhile,  passive  optical  limiters  rely  solely  on  the  incoming  laser 
radiation  for  switch  activation.  As  a  result,  the  switching  time  can  be  minimized  and 
very  little,  if  any,  auxiliary  equipment  is  required. 

Passive  optical  limiters  can  be  constructed  by  using  a  nonlinear  optical  material 
in  a  specific  geometrical  configuration.  A  nonlinear  material  is  one  which  displays  an 
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index  of  refraction  or  an  absorption  which  is  dependent  on  the  incident  laser  intensity. 
The  deflection  of  radiation  or  the  absorption  of  radiation  can  thus  increase  as  the 
incident  laser  radiation  increases.  This  increase  in  energy  absorbed  or  deflected 
produces  a  decrease  in  transmission  and  thereby  provides  self-generated  protection 
against  a  laser  threat. 

The  effort  described  here  will  concentrate  on  a  dynamic  passive  device  via 
photorefractive  beam  fanning.  The  device  operates  on  the  principle  of  dynamic 
gratings.  In  the  formation  of  a  dynamic  grating  the  incoming  laser  radiation  produces  a 
modulation  in  the  material  index  of  refraction.  The  index  modulation  then  acts  as  a  real 
time  phase  grating  and  deflects,  or  scatters,  the  incident  laser  light  as  depicted  in  Figure 
1-1. 

There  are  many  benefits  which  can  be  derived  from  the  development  of 
photorefractive  dynamic  grating  devices.  For  example,  they  offer  the  following 
desirable  features  (Figure  1-2): 

(i)  low  switching  threshold  to  protect  against  jamming  and 
spoofing; 

(ii)  high  optical  density  (OD)  to  protect  against  laser-induced 
damage; 

( Hi)  broad  spectral  bandwidth  to  cover  in-band  agile  threats; 

( iv)  can  protect  against  short  plusses  well  as  cw  radiation ; 

(v)  direction  insensitivity  to  protect  against  off-angle  radiation; 

(vi)  allow  the  outside  world  to  be  monitored  while  it 
simultaneously  rejects  laser  radiation. 

Lets  look  at  one  of  these  features  more  closely.  A  photorefractive  limiter  has 
the  potential  to  protect  against  both  short  pulse  and  c.w.  threats.  This  capability  is  a 
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SIMPLE  LIMITER 


result  of  the  fact  that  the  laser  light  intensity  transmitted  through  a  photorefractive 
limiter,  IT  falls  off  nearly  exponentially  in  time,  i.e. 


IT=I0e-,/T  (1-1) 

where  the  characteristic  response  time  is  found  to  be  inversely  proportional  to  the 
incident  intensity,  IQ,  and  given  by 


X 


(1-2) 


A  is  a  constant. 

As  a  result,  the  energy  per  unit  area  which  passes  through  the  crystal  is 
therefore  given  as 


J  ITdt  =  Jl0e"t/X  dt 


d-3) 


Using  (1-2)  in  (1-3)  we  find 

~  T 

[  ITdt=— -=A.  (1-4) 

J  X 

o  * 


Experimental  values  of  the  constant  “A”  are  on  the  order  of .  1  to  .01  Joules  per  cm2.  It 
is  immediately  obvious  that  if  “A”  could  be  reduced  by  about  four  to  five  orders  of 
magnitude,  the  energy  transmitted  trough  the  beam  fanning  limiter  would  be  clamped  in 
the  microjoule  range.  Amazingly  this  limiting  value  would  be,  therefore,  independent 
of  the  input  temporal  width  or  input  intensity.  Physically,  this  attractive  feature  is  due 
to  the  fact  that  if  the  incident  intensity  increases  the  photorefractive  limiter  response  time 
decreases  proportionally,  switching  off  the  higher  incident  intensity  more  quickly.  The 
faster  shutter  speed  results  in  the  transmitted  energy  remaining  constant  at  the  value 
“A”.  Of  course,  in  practice,  the  transmitted  energy  does  not  drop  to  zero  but  rather  to 
an  equilibrium  value  I  .  The  comparison  between  ID  and  I  is  called  the  optical 
density,  or  O.D.,  of  the  limiter.  That  is. 
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(1-5) 


O-D.slog— 

U 

The  steady  O.D.  is  determined  by  the  photorefractive  index  change  which 
depends  on  the  magnitude  of  the  effective  electro-optic  coefficient  of  the  crystal. 
Presently  available  materials,  such  as  BaTi03  and  SBN,  have  large  electro-optic 
coefficients  and  consequently  display  strong  photorefractive  index  changes.  When  this 
strong  photorefractive  index  change  is  maintained  by  preserving  the  gain-length 
product,  then  significant  improvements  in  the  performance  of  beam  fanning  limiters  can 
be  realized  by  using  a  number  of  different  techniques  to  enhance  (speed-up)  the  decay 
time  of  the  transmitted  beam. 

Any  technique  employed  to  improve  the  beam  fanning  response  time,  tBF,  must 
not  be  accompanied  by  a  loss  of  gain  or  the  value  of  the  O.D.  will  suffer.  Optical 
densities  in  the  steady-state  should  approach  3  or  4.  Although  the  parameter  A  has  not 
been  studied  in  great  detail,  however,  it  is  known  to  be  a  function  of  dopants, 
wavelength,  temperature,  focusing  geometry,  length,  seeding,  applied  electric  field, 
and  surface  quality.  A  number  of  techniques  have  been  devised  recently  for  improving 
the  response  time  of  photorefractive  materials  without  sacrificing  gain.  These  include 
for  example,  the  techniques  of  increasing  the  laser  intensity,  using  applied  electric 
fields,  seeding  the  scattering,  causing  self-focusing  and  defocusing,  controlling  the 
dopant  concentration  and  changing  the  crystal  temperature.  To  discuss  the  merit  of  the 
results  of  our  investigation  of  these  enhancement  techniques  we  now  discuss  a  model 
for  the  photorefraction  limiter. 

Conclusion 

It  is  clear  that  our  goal  is  to  improve  the  photorefractive  response  time  by  4  to  5 
orders  of  magnitude  while  maintaining  an  O.D.  or  3  to  4. 


7 


II. 


PHOTOREFRACTION 


The  most  widely  accepted  model  used  to  explain  the  photorefractive  effect  is  the 
band  transport  model.  This  model  assumes  the  existence  of  mobile  charge  carriers 
(which  are  usually  electrons  in  SBN)  residing  in  traps  in  the  material.  When  an 
intensity  distribution  is  incident  on  the  material,  the  charge  carriers  in  the  bright  regions 
of  the  material  are  photoexcited  from  the  traps  into  the  conduction  band.  These  excited 
charges  then  migrate  via  drift  and  diffusion  leaving  behind  ionized  traps.  The  excited 
charge  carriers  are  continually  retrapped  and  excited  again  as  they  migrate  until  they  are 
eventually  retrapped  in  thr  dark  regions  of  the  material,  where  the  rate  of 
photoexcitation  is  greatly  reduced.  As  a  result  of  this  migration,  a  charge  distribution 
forms  in  the  material  consisting  of  the  ionized  traps  in  the  bright  regions  and  the 
retrapped  charge  carriers  in  the  dark  regions.  This  charge  distribution  generates  a 
space-charge  field  in  the  material  in  accordance  with  Gauss’s  law.  The  space-charge 
field,  in  turn,  produces  a  change  in  the  refractive  index  of  the  material  through  the  linear 
electro-optic  (PockePs)  effect. 

If  two  plane  waves  interfere  coherently  in  the  material  so  that  the  intensity 
distribution  varies  periodically,  then  the  resulting  index  change  will  be  periodic  as  well. 
Thus,  an  index  grating  is  produced  in  the  material.  According  to  Gauss’s  law,  there 
will  be  a  jt/2  phase  shift  between  the  incident  intensity  distribution  and  the  induced 
index  grating.  The  tc/2  phase  shift  will  cause  the  two  waves  to  couple  as  the  propagate 
through  the  material  and  will  transfer  energy  from  one  wave  to  the  other  in  the  process. 
In  general,  this  phase  shift  between  the  incident  intensity  distribution  and  the  induced 
index  grating  is  uniquely  responsible  for  phenomena  which  make  photorefractive 
materials  inviting  as  an  optical  limiter. 

One  of  the  goals  of  the  mathematical  analysis  of  the  band  transport  model  is  an 
expression  for  the  space-charge  field  in  the  photorefractive  material.  A  preliminary 
framework  must  be  established  in  order  to  analyze  the  problem.  The  positive  x-axis 
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Fig.  3-4  The  steps  composing  the  photorefractive  effect 


will  be  defined  by  the  c-axis  of  the  photorefractive  material.  Two  monochromatic  plane 
waves  of  the  same  frequency  propagate  through  the  material  with  wave  vectors  it,  and 
k2 .  These  plane  waves  interfere  coherently  to  produce  a  sinusoidal  intensity  pattern 
whose  grating  wavevector  is  kg  =  kt  -  k2.  To  simplify  the  problem,  the  following 
constraints  are  placed  on  the  wavevectors  of  the  plane  waves:  KyI  =  ky2  =  0,  kzl> 0  and 
kzp> 0,  i.e.  the  propagate  in  the  xz  plane. 

Two  species  of  impurity  atoms  are  present  in  the  photorefractive  material.  The 
density  of  acceptor  impurities  is  given  by  NA  and  is  assumed  constant.  The  density  of 
donor  impurities  is  given  by  ND  likewise  considered  constant.  The  energy  state  of  the 
acceptors  is  assumed  to  be  lower  than  that  of  the  donors.  Thus,  each  acceptor  takes  an 
electron  from  a  donor,  producing  an  ionized  donor  capable  of  trapping  excited  charge 
carriers.  It  is  assumed  that  the  acceptor  impurities  NA  are  fully  occupied  at  all  times,  so 
they  do  not  participate  actively  in  the  charge  transport  mechanism.  The  density  of 
ionized  donors  is  denoted  by  N+D(x,t\  while  the  density  of  electrons  in  the  conduction 
band  is  given  by  ne  (x,t).  When  no  light  is  incident  on  the  material,  njx,t)  =  0  and 

R)=na. 
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Rate  Equation  :  the  rate  of  change  in  the  number  of  ionized  donor  atoms  (neglecting 
thermal  ionization)  is  the  difference  in  the  rate  at  which  the  donors  are  excited  and  the 
rate  at  which  ionized  donors  recapture  electrons: 

3A It 

Jt  =  si  (x)  (Nd  -  N*d)  -  yRne ,  (2.1) 

where  5  is  the  photoionization  cross-section,  I(x)  is  the  incident  intensity  distribution 
along  the  x-axis  and  yR  is  the  recombination  coefficient  of  a  free  electron  at  the  site  of 
an  ionized  donor. 

Current  Density  :  the  (conventional)  current  density,  jx  ( x ,  t ) ,  is  composed  of  drift  and 
diffusion  terms: 


9 n 

/x  ~  Q en > 


(2.2) 


where  qe  is  the  magnitude  of  electron  charge,  \xe  is  the  electron  mobility,  Ex  {x,  t)  is  the 
total  electric  field  (internal  and  external),  kB  is  Boltzmann’s  constant  and  T  is  the  tem¬ 
perature  (in  Kelvin). 


Continuity  Equation  :  any  spatial  variation  of  the  current  density  must  be  accompanied 
by  the  accumulation  or  depletion  of  charge  carriers  in  the  region  of  variation: 


(2.3) 


Gauss’s  Law  :  the  spatial  variation  of  the  electric  field  is  dependent  upon  the  charge  dis¬ 
tribution  in  the  material: 


dE 

X 

dx 


(N*D-nt-NA)  , 


where  e  is  the  static  permittivity  of  the  material. 
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(2.4) 


Solution  for  the  Space-Charge  Field 

The  equations  in  the  preceding  section  form  a  highly  nonlinear  system  which,  in  general, 
cannot  be  solved  analytically.  However,  a  closed  form  solution  can  be  obtained  with  some 
assumptions.  The  interference  of  two  monochromatic  plane  waves  results  in  an  intensity 
distribution  which  is  modulated  sinusoidally.  Hence,  if  I(x)  is  expressed  as  a  Fourier  series 
and  it  is  assumed  that  all  but  the  zeroth-  and  first-order  terms  are  negligible,  then  the  inten¬ 
sity  distribution  can  be  written  as  follows: 

/(*)  =  IQ+^UleikgX  +  c.c.]  =  I0+j[I0eik*x  +  c.c.]  ,  (2.5*) 

where  m  =  Zj/Zq is  the  modulation  ratio.  The  other  variables  then  assume  similar  forms: 


ne  (x,  t)  -  neQ(t)  +^[nel(t)  e‘kgX  +  c.c.]  ’  (2.5.b) 

A£(x,  0  =  (0  +  i  [ATf  (/)  e‘kgX+c.c.)  ,  (2.5.c) 

ix{x,  t)  =  jQ (0  + 1  Ui  (0 eikgX  +  c.c.]  ,  (2.5. d) 

Ex(x,t)  =  E0(t)  [El  (t)  e'kgX  +  c.c.]  ,  (2.5.e) 


where  Eq  is  any  applied  external  field  and  Ex  is  the  space-charge  field  in  the  photorefrac- 
tive  material.  Substituting  these  equations  into  (2.1)-(2.4),  the  following  set  of  linearized 
equations  is  obtained: 

dNt 

Jt  =  sIo  Wd-ND-  yRne0N^  ,  (2.6.a) 


dNt 

Jt  =  smIQ(ND-N+) -slQl^-y^n^  +  n^) 
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Jo  =  <leVene  0E0 


(2.6.C) 


h  =  'We(ne0El+nelE0 )  ’ 


0  = 


<K_dneQ 

dt  dt 


'Vi  =  -^U  -*  J- 


0  -  Hq  -  neo~ Na 


'Vl  =  f  K""el)  • 


With  these,  one  can  obtain  the  following  second-order  differential  equation  for  E\ : 


d~E 


1  .  ,  1  .1  .  1  i  dE\  .  .  1  1 


,  1  / 


3J  +(r+;r+rT)^  +(;rr+;rr_?T)£i  =  7Fi(r~r)™,/o(A'D~n‘o_A^) 

dt  l+  ldi  '■D  E  +  Vi  TD  X/X£  'SE  XD  XE 


where  the  time  constants  in  (2.7)  are 


q 

xn  =  mean  diffusion  time  =  - -  , 

\i  kRTk\ 

r e  B  g 


Trs  mean  drift  time  =  — , 

M»£  o 


%di  =  dielectric  relaxation  time  = 


x.  =  excitation  time  =  — — i -  , 

1  sI0  +  yRne0 


(2.6.d) 

(2.6.e) 

(2.6.f) 

(2.6.g) 

(2.6.h) 

,  (2.7) 

(2.8.a) 

(2.8.b) 

(2.8.c) 

(2.8.d) 
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1 


(2.8.e) 


x.  = 


+  sI0  +  2yRne0  +  yRNA 


For  low  irradiances,  i.e.  few  photoexcited  carriers  (NA » ng0)  and  fast  trapping 
(jRNA  »  sIq),  x+  reduces  to  xR,  the  recombination  time: 


x.  ~ 


*  1rNA 


=  X 


R  ■ 


(2.9) 


A  detailed  derivation  of  (2.7)  is  given  in  Appendix  I,  as  well  as  derivations  of  the  steady- 
state  and  time-dependent  forms  of  the  space-charge  field  in  the  photorefractive  material. 

Steady-State  Solution 

The  steady-state  solution  for  the  space-charge  field  is  easy  to  obtain  since  the  time  deriva¬ 
tives  in  (2.7)  vanish.  For  low  irradiances,  the  equilibrium  space-charge  field  is  given  by 


E?  =  mEj* , 


(2.10) 


where  the  amplitude,  Esc,  and  the  phase  shift,  4>,  are 


Esc  ~  en. 


E2d  +  E o 


(£0 +  £„)*  +  £* 


(2.11. a) 


f 


tan  ({)  = 


1  +  ^  +  ^ 
EN  edeN/ 


(2.1  l.b) 


and  Ed  and  E ^  are  defined  as 


kRTk0 

p  -  _1 _ | 


(2.12.a) 


E  ~  qeNgff 

N~  V 


1eNA  ( 
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Na  \ 


k+»d-"aj 


-1 


(2.12.b) 


The  electric  field  Eq  is  the  effective  field  causing  electron  diffusion  while  the  electric  field 
En  is  the  maximum  space-charge  field  attainable  without  applying  an  external  field.  The 
total  electric  field  at  equilibrium  is  then 


EsJ(x)  =  E0  +  ^Esc[ei(kgX+^  +  c.c.]  . 


(2.13) 


Time-Dependent  Solution 

Due  to  the  form  of  (2.7),  some  damped,  oscillatory  behavior  in  the  space-charge  field 
would  be  expected  as  it  evolves  from  its  initial  value  to  its  equilibrium  value.  If  no  gratings 
exist  in  the  photorefractive  material  at  time  t  =  0  and  if  we  assume  that  x+  «  xdj,  x+  «  %D 
and  x+  «  Xj,  then  the  time-dependent  form  of  the  space-charge  field  is  given  by 


- t/z  /cor-, 


Ex{t)  =  E\s[\-e  "V1"']  , 


(2.14) 


where  x  and  ©  are 


1+Ii  +  M2+f!if 

XD  Xdi)  VX£y 


XT  .. 
+  di 


i  +  —  +  —+  —  [  — 

XD  Xdi)  XAXE 


(2.15.a) 


xA(xdi 


1+!i+i±f+r!iV 

XD  Xdi  7  \XE  / 


(2.15.b) 


The  time-dependent  form  of  the  total  electric  field  is  then 


Ex(x,t)  =  E0+^Esc[e‘^x+^  [\-e',/xeim]  +c.c.] 


(2.16) 
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Pockel’s  Effect  and  Grating  Formation 

The  space-charge  field  produced  by  charge  separation  in  the  photorefractive  crystal 
induces  a  change  in  the  refractive  index  of  the  crystal  via  Pockel’s  effect.  Given  an  electric 
field  in  the  material,  the  change  in  the  impermeability  tensor  is  given  by 

Al1  U  =  rijkEk  -  (2.17) 

where  r ^  is  the  linear  electro-optic  tensor.  This  results  in  an  index  change  given  by 

An  =  ~\nlreffE  ,  (2.18) 

where  n0  is  the  unperturbed  refractive  index  of  the  medium  and  rejj  is  a  linear  combination 
of  elements  of  r ^  that  depends  on  both  orientation  of  the  crystal  and  direction  of  the  elec¬ 
tric  field  (see  Appendix  II).  By  substituting  (2.13)  into  (2.18),  we  obtain  the  following 
expression  for  the  equilibrium  change  in  the  refractive  index  due  to  the  electric  field  in  the 
material: 


1  3 

—  _  _ n  i 


An  =  -An^ffE0-^nlrlffEsc  +  c.c.  ]  . 


/(*-*+$) 


(2.19) 


The  spatial  variation  of  the  refractive  index  in  the  material  can  be  expressed  as 


/  \  1  r  IJtjT  , 

n(x)  =  nl  +  -[n2e  *  +c.c.]  . 


(2.20) 


By  comparing  (2.19)  and  (2.20),  n\  and  «2  must  be  given  by  the  following: 


ni  =  n0  ~  \nlreffE0  ’ 


(2.21. a) 


m  3  -  id 

n2  =  ~2n0  reffEsce 


(2.2  l.b) 
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Assuming  reff  is  positive,  (2.21.b)  shows  that  the  modulation  of  the  refractive  index  is 
in  phase  with  the  space-charge  field  and,  thus,  shifted  by  the  phase  angle  <f>  with 
respect  to  the  incident  intensity  distribution. 

Summary 

Having  coumpleted  a  description  of  the  basic  underlying  physics  for  the 
photorefractive  limiter  we  can  now  turn  our  attention  to  a  discussion  of  the  most 
significant  results  of  our  efforts  to  enhance  the  photorefractive  response  time  without  a 
reduction  in  the  O.D. 
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III.  INCREASING  THE  INTENSITY 


The  photorefractive  response  time  is  given  in  expression  (2.15.a).  For  our 
materials,  x+  ~  xR  «  x^x^  ,xt .  a  result,  for  zero  applied  Field  (2. 15.a)  becomes 
x=xdi  (3-1) 

Since  xd  is  inversely  proportional  to  neo  which  is  itself  proportional  to  the  incident  total 
intensity  I0,  x*  is  inversely  proportional  to  IQ,  i.e. 


Physically,  this  is  the  case  because  a  higher  intensity  means  more  charges  in  the 
conduction  band  and  therefore  an  increased  diffusion  rate  and  a  faster  approach  to 
equilibrium. 

Given  that  the  photorefractive  response  time  can  increase  due  to  a  higher 
intensity  it  is  easy  to  see  that  one  approach  to  achieving  our  goal  of  a  faster  response 
time  would  be  to  focus  the  incident  laser  power  to  a  small  spot  size.  In  this  way,  for  a 
given  incident  power  the  photorefractive  response  time  will  shorten  dramatically.  This 
approach,  however,  is  not  straight  forward. 

The  beam  fan  originates  from  conventional  scattered  light  which  crosses  the 

incident  beam  and  writes  index  gratings  in  the  crystal  via  the  photorefractive  effect.  As 

a  result  of  these  index  gratings,  energy  is  coupled  from  the  incident  beam  into  the 

weaker  scattered  beams,  which  can  result  in  significant  depletion  of  the  incident  beam. 

The  magnitude  of  the  coupling  depends  on  the  number  of  scatterers,  the  gain  of  the 

crystal,  and  the  available  length  ((interaction  length)  over  which  the  incident  beam  and 

the  scattered  beam  the  interaction  length  can  typically  be  several  millimeters  long  (i.e., 

of  the  order  of  a  few  beam  waists).  This  interaction  length  does  not  change  as  the 

intensity  is  increased  by  turning  up  the  laser  power  and  an  intensity  dependence  of 

approximately  I'1  is  typically  observed  for  the  beam  fanning  time  response.  However, 

for  abeam  focused  with  a  spherical  lens  the  interaction  length  may  be  of  the  order  50 
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1  W/cm2 


INTERACTION  LENGTH 


Interaction  length  or  Interaction  length  or 

7L  mostly  large  yL  mostly  small 


(im  or  less.  Beam  fanning  is  almost  absent  for  such  small  interaction  lengths  even 
though  the  time  response  still  exhibits  an  I'1  intensity  dependence. 

The  situation  is  different  with  a  cylindrical  lens  which  can  produce  faster 
response  times  without  a  reduction  in  coupling  strength.  This  is  because  a  cylindrical 
lens  focuses  in  one  dimension.  It  can,  therefore,  produce  a  tight  focus  in  the  dimension 
which  has  a  small  gain  and  still  maintain  a  large  interaction  length  in  the  other 
dimension.  For  example,  in  SBN  the  relevant  EO  coefficients  are  r33  =  2224  x  10' 
l2m/V  and  r13  =  55  x  10'12m/V,  where  r33  provides  the  maximum  gain  for  extraordinary 
polarized  beams  propagating  at  near  normal  incidence  to  the  c-axis  (the  grating  wave 
vector  parallel  with  the  c-axis).  As  a  result,  the  greatest  beam  depletion  is  achieved  by 
combining  the  large  gain  associated  with  r33  and  the  large  interaction  length  associated 
with  the  major  axis  of  the  elliptical  beam  spot  produced  by  the  cylindrical  lens.  The 
gain  coefficient  in  the  other  direction  is  dominated  by  rI3  and  is  small  by  comparison, 
so  little  is  lost  by  having  a  tight  focus  in  this  direction.  In  the  cylindrical  focusing 
geometry  the  intensity  is  increased  although  not  as  much  as  the  spherical  lens  case  (i.e., 
Jt  i2,,  vs  itr0r).  However,  the  increase  in  intensity  decreases  the  response  time  at  no  cost 
to  beam  depletion. 

A  comparison  of  the  beam  fanning  response  is  shown  in  Fig.  3.1.  In  this 
figure  the  c-axis  of  the  crystal  is  into  the  plane  of  the  page  and  the  beams  are  focused  in 
the  crystal.  In  these  conditions  there  is  virtually  no  beam  fanning  produced  when 
focusing  to  a  spot  size  of  0. 1mm  (FWl/e2M)  with  a  spherical  lens.  This  comparison  is 
somewhat  striking  since  the  unfocused  spot  size  was  2.12mm,  and  the  cylindrical  lens 
focused  the  beam  to  an  elliptical  spot  with  axes  of  0. 1 1  x  2. 12mm. 

The  response  times  vs  power  for  unfocused  and  cylindrical  focusing  are 

compared  in  Fig.  3.2  for  the  same  spot  sizes  used  in  the  measurements  of  Fig.  3.1. 

The  observed  shortening  of  the  response  time  (factor  of  20)  by  the  cylindrical  lens 

compares  well  with  the  ratio  of  the  areas  of  the  focused  and  unfocused  beams.  In  other 
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words,  the  improvement  in  the  time  response  obtained  for  the  cylindrical  lens  is  just 
due  to  the  increase  in  intensity  achieved  by  reducing  the  spot  size.  The  parameters  A 
and  b,  given  in  the  caption,  were  obtained  from  a  least-squares  fit  to  the  data  and  agree 
well  with  previous  measurements  of  the  beam  fanning  response  time  where 
approximate  I"1  intensity'  dependence  was  found  for  Ce  doped  SBN  crystals.  The 
present  measurements  indicate  that  the  I' 1  intensity  dependence  of  the  time  response  is 
valid  over  at  least  3  orders  of  intensity  variation. 

Conclusion 

We  have  demonstrated  that  the  use  of  a  cylindrical  lens  to  focus  the  incident 
laser  power  down  to  a  small  spot  size  can  reduce  the  response  time  while  maintaining 
an  O.D.  of  2  to  3.  This  means  that  an  incident  beam  of  5cm  x  5cm  can  be  focused 
using  a  cylindrical  lens  50  microns  x  5cm  or  a  reduction  in  size,  in  one  dimension,  of 
10‘3.  This  results  in  an  intensity  increase  of  103  and  a  corresponding  decrease  in 
response  time  of  3  orders  of  magnitude.  Since  this  is  close  to  our  4  to  5  orders  of 
magnitude  goal,  it  is  clear  that  this  technique  is  significant  when  coupled  with  others 
may  produce  a  viable  photorefractive  limiter. 
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TIGHT  FOCUS 


BEAM  FANNING  AT  THE  FOCAL  POINT 
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TIME  (sec) 


BEAM  FANNING  RESPONSE  TIME 


IV.  APPLIED  DC  ELECTRIC  FIELD 


As  can  be  seen  from  expressions  (2.1  l.a)  and  (2.15.a),  both  the  magnitude  of 
the  space  charge  electric  field  and  its  time  of  formation  response  time  can  be  enhanced 
with  the  application  of  a  d.c.  electric  field. 

The  dependence  of  the  space  charge  field  and  its  formation  time  on  the  applied 
field  can  be  seen  from 


approaching  =  -90°  for  large  applied  fields 

and  T  ->  I,  (4-5) 

approaching  x,  for  large  applied  fields 
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and 
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(4-6) 
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approaching  a  number  which  make  possible  to  observation  of  oscillations. 

Physically,  Eq  or  is  the  electric  field  that  would  exist  if  all  available  charge 
carrier  bright  to  the  dark  regions  of  an  interference  pattern.  Meanwhile,  x,  is  the  time 
to  excite  ah  the  charges  will  be  transported.  In  addition,  the  period  of  oscillation  T  is 
now  the  ratio  of  xE  to  xR  times  the  response  time  x,.  In  other  words  it  is  related  to  the 
number  of  grating  periods  the  carriers  travel  before  re-trapped. 

While  these  limiting  values  provide  a  sense  of  the  potential  of  using  an  applied 
electric  field,  numerical  solutions  can  be  used  to  focus  more  at  finding  optimum  crystal 
parameters. 

The  numerical  method  used  in  our  computer  program  and  the  physical 
parameters  are  listed  in  Appendix  IV  and  V.  The  computer  program  was  run  for 
different  trap  number  densities  of  the  crystal.  For  each  value  of  NA  the  total  space- 
charge  field  was  calculated  for  both  the  no  field  case  and  the  DC  field  case.  The 
calculated  results  for  the  space-charge  field  time  evolution  are  shown  in  Fig.  4-1 
through  Fig.  4-6. 
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These  curves  show  that  due  to  a  DC  field  (1)  the  magnitude  of  the  space-charge  field 
was  increased;  (2)  the  photorefractive  response  time  was  reduced;  and  (3)  an 
oscillation  of  the  space-charge  field  appeared  ( see  Fig.4-6 ). 

For  the  curves  Fig.4-1  through  Fig.4-6  all  the  parameters  are  the  same  except 
the  trap  number  density  in  the  crystal.  The  NA  values  picked  in  different  curves  are 
5*1016cm’3  in  Fig.4-1  and  4-2,  1.5*  1016  cm’3  in  Fig.4-3  and  4-4,  and  5*1015  cm'3  in 
Fig.4-5  and  4-6.  By  analyzing  and  comparing  these  curves  we  found  that  for  crystals 
with  different  NA  values,  the  enhancement  of  the  photorefractive  effect  was  different 
even  with  the  same  DC  field.  Looking  at  Fig.4-1  and  Fig.4-2,  it  is  seen  that  for  a  large 
Na  value  the  applied  DC  field  doesn’t  change  the  photorefractive  response  time 
substantially.  However,  by  reducing  the  NA  value  from  5*1016cm’3to  1.5«1016cm'3  , 
the  results  shown  in  Fig.4-3  and  Fig.4-4  are  exciting.  We  see  that  with  the  same  DC 
field  as  used  in  the  above  case  a  larger  enhancement  is  observed  in  both  the  magnitude 
and  the  speed  of  the  photorefractive  effect.  The  response  time  is  reduced  by  about  one 
half  with  a  DC  field  of  4  kv/cm  to  a  crystal.  In  Fig.4-5  and  Fig.4-6  the  NA  value  was 
decreased  further  to  5*1015  cm’3.  We  see  that  the  response  time  is  also  reduced  to 
about  one  half  with  a  DC  field  of  4  kv/cm,  but  the  actual  response  time  in  this  case  is 
longer  than  that  in  the  case  of  NA=1.5*1016cm3.  Moreover,  an  oscillation  with  time  in 
the  space-charge  field  is  observed.  From  eq.(4-6)  we  can  get  a  dependence  of  co  on  NA 
and  Nd  ,  which  is 

£0  a  (  ND2  -  Na2)/[  c,(  Na2+  c2 )  +  c3  ]  ( 4-7  ) 
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where  Cj,  c2  and  c3  are  constants  independent  of  NA  and  ND.  From  (4-7)  we  can  see 
that  a  smaller  NA  will  cause  a  bigger  co.  Physically,  a  smaller  NA  makes  a  longer 
recombination  time  xR,  which  results  in  that  the  electrons  in  conduction  band  can  move 
freely  across  one  or  more  grating  space  periods  without  being  recombined.  That’s  why 
Fig.4-6  showed  an  oscillation. 

Therefore,  the  conclusion  is  that  if  the  ND  value  is  fixed  an  appropriate  NA  value 
of  a  crystal  is  vital  to  observe  an  enhancement  of  the  photorefractive  response  time 
with  an  applied  field.  This  important  result  can  be  explained  physically.  The 
recombination  time  xR  is  inversely  proportional  to  the  crystal  trap  number  density. 
With  Na  decreasing  xR  is  increasing.  As  xR  is  increased  to  be  longer  than  the  field  drift 
time  Xg,  electrons  in  the  conduction  band  can  move  further  without  being  recombined. 
A  longer  recombination  time  makes  it  possible  for  the  space-charge  field  to  be  built  up 
in  a  shorter  time.  This  is  why  the  photorefractive  response  time  is  getting  shorter  for 
the  same  magnitude  of  an  applied  field  as  the  crystal  has  a  smaller  NA  value.  However, 
if  Na  is  too  small,  or  in  other  words,  the  recombination  time  xR  is  too  long,  electrons  in 
conduction  band  can  move  freely  across  several  grating  space  periods  without 
recombination.  This  will  result  in  an  oscillation  of  the  space-charge  field,  and  it  will 
take  longer  time  to  have  the  space-charge  field  reach  an  equilibrium.  So  it  seems  that 
there  is  a  compromise  for  a  trap  number  density  of  the  crystal  between  a  big  number 
and  a  small  number.  The  value  of  NA  should  be  chosen  such  that  by  applying  an 
external  electric  field  across  a  crystal  the  resulted  space-charge  field  will  be  able  to 
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reach  an  equilibrium  faster  but  without  causing  oscillations.  We  also  ran  the  computer 
program  for  the  space-charge  field  by  keeping  NA  constant  while  changing  the  ND 
value  and  kg  value,  respectively.  The  parameters  and  results  are  as  follows: 

1.  Varying  ND  value 


DC  field:  4  kv/cm  NA=  1.5*1016cnr3  ND=  5*1017  cnr3 
kg  =  2.4‘IQ4 1/cm  I  =  100  mw/cm2 


Fig.  4-7 

DC  field:  4  kv/cm  NA=  1.5*1 01 6 cnr3  ND=  1.5*1017cm'3 
(smaller  Nd)  kg  =  2.4*1 04 1 /cm  I  =100  mw/cm2 


Fig.  4-8 
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DC  field:  4  kv/cm  NA=  1.5*1016cm'3  ND=  1.5*10'8cnr3 
(bigger  ND)  kg  =  2.4*1 04  1/cm  I=100mw/cm2 


The  curves  showed  that  (1 )  a  smaller  ND  made  a  longer  response  time;  (2)  a 
larger  N0 caused  an  oscillation  of  the  space-charge  field.  Mathematically,  from 
expression  (4-7)  we  see  that  a  bigger  ND  causes  a  bigger  co.  Physically,  with 
same  optical  intensity,  if  NDis  bigger,  more  electrons  will  be  excited  to  the 
conduction  band  than  the  case  with  smaller  ND.  This  will  make  more  electrons  in 
conduction  band  moving  freely  across  one  or  more  grating  space  periods. 
Therefore,  an  oscillatory  space-charge  field  is  produced. 

2.  Varying  k„  value 

In  Fig.4-7,  Na=1.5*1016  cm'3,  ND=5*1017  cm*3,  kg=2A*\W  cm-'; 

In  Fig.4-10,  Na=1.5*1016  cm'3,  ND=5*10n  cm*3,  kg=8*103  cm'1  ( smaller  kg ); 

In  Fig.4-11,  Na=1.5*1016  cm'3,  ND=5*1017  cm-3,  kg=7.2*104  cm'1  ( bigger  kg ). 
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DC  field:  4  kv/cm 


NA  =  1.5*1016cm-3  ND  =  5*1017cm-3 


smaller  k  (k=8.0*  10  ^  1/cm) 


Fig.  4-10 


DC  field:  4  kv/cm 


NA  =  1.5‘IQ1^  cm'3  ND  =  5*1017  cm'3 


200  400  600  800  1000  t  (xO.Ol  s) 

Fig.  4-11 


The  curves  showed  that  (1)  a  smaller  kg  produced  a  longer  response  time  ;  (2)  a 
bigger  kg  produced  a  shorter  response  time  but  a  smaller  enhancement  in  the 
magnitude  of  the  space-charge  field.  We  plotted  the  curves  of  x  versus  ND  and 
versus  kg,  which  are  shown  in  Fig.4-12  and  Fig.4-13. 


u  5  10  15  20  25  30 

Nd (xlO17  cm-3) 

Fig.  4-12  Dependence  of  Response  Time  on  ND 


x(s)  NA=  1.5«1016cnr3  ND=5*10]7  cm-3 


2  4  6  8  10 

kg  (xlO4  1/cm) 

Fig.  4-13  Dependence  of  Response  Time  on  kg 


From  Fig.4-12  we  can  see  that  a  bigger  ND  makes  a  shorter  t,  and  from  Fig.4-13 
we  see  that  a  bigger  kg  makes  a  shorter  x.  The  physical  explanation  is  that  (1)  A 
bigger  ND  makes  a  shorter  x,,  thereby,  produces  a  shorter  response  time;  (2)  A 
bigger  kg  makes  a  shorter  x£.  As  xE  gets  approximately  equal  to  or  less  than  xR, 
the  response  time  will  be  reduced.  So  from  these  two  curves  we  are  able  to 
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predict  the  optimum  ND  and  k„  values.  For  a  DC  field  of  4  kv/cm,  if  NA  value 
was  more  or  less  than  1.5*1016  cm'3,  the  numerical  results  showed  that  the 
response  time  was  longer  than  the  case  with  NA  value  of  1.5M016  cm'3. 
Therefore,  for  a  DC  field  of  4  kv/cm  and  ^=10’9cm3/s,  the  NA  value  of  1.5*1016 
cm'3  is  the  optimum  value  for  getting  a  fast  response  time,  and  the  optimum  ND 
and  kg  values  are  5*1017cm’3and  2.4*104cm’1,  respectively.  In  the  general  case, 
for  a  DC  field  of  4  kv/cm  the  optimum  NA  value  is  about  1 .5*  1 07/yR  cm'3,  and  the 
curve  of  response  time  versus  ND  and  the  curve  of  response  time  versus  kg  are 
similar  to  Fig.4-12  and  Fig.4-13.  Therefore,  we  can  estimate  that  the  optimum 
ND value  is  about  30  times  of  NA  value  and  the  optimum  kg  value  is  in  the  range 
of  2*104  cm'1  to  3«104  cm*1.  In  other  words,  the  optimum  full  crossing  angle 
between  two  beams  is  in  the  range  of  8°  to  13°. 

Experimental  Setup  And  Data  Presentation 

Experimentally,  we  measured  the  photorefractive  signals  by  applying  a  DC  field 
across  a  crystal.  Our  experimental  research  can  be  described  as  follows. 

(1) .  The  measurement  technique  adopted  in  the  experiment  included  beam  fanning 

and  two  beam  coupling. 

(2) .  Samples  measured  by  us  included  a  heavily  doped  BSKNN  crystal  which  had  a 

larger  value  of  NDand  NA,  and  an  undoped  BSKNN  crystal  which  had  a  smaller 
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value  of  Nd  and  NA 

(3).  Each  sample  was  measured  with  no  field  and  with  an  applied  field  across  it . 


Experimental  Apparatus 


The  experimental  setup  for  beam  fanning  is  shown  in  Fig.4-14.  The  laser  we  used 


Digitizing 


Polorizer 


Fig.  4-14  Experimental  setup  for  beam  fanning 

experiment  was  a  Coherent  Argon-ion  Laser.  The  laser  beam  of  458  nm  wavelength 
with  an  extraordinary  polarization  was  incident  to  a  crystal  in  the  direction 
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perpendicular  to  the  c-axis  of  the  crystal.  The  transmitted  beam  was  detected  by  a 
photodiode  through  an  aperture.  The  signal  was  then  sent  to  a  602A  Digitizing  Signal 
Analyzer  to  be  measured.  In  order  to  measure  the  photorefractive  effect  with  a  DC 
field  a  high  voltage  needs  to  be  applied  across  a  crystal.  Two  copper  tapes  were  pasted 
on  a  crystal  by  using  silver  paste  and  the  electrodes  were  connected  to  a  high  voltage 


supply.  A  high  voltage  was  provided  by  a  TREK  Model  610C  H.V.  SUPPLY 
AMPLIFIER/CONTROLLER  which  had  a  maximum  output  of  10  kv/cm. 


The  setup  for  two  beam  coupling  is  shown  in  Fig.4-15.  The  laser  output  with 


ordinary  polarization  and  a  wavelength  of  458  nm  was  split  into  two  beams.  By  using 
neutral  density  filters  we  made  the  intensity  of  one  beam  much  stronger  than  the  other. 
The  full  crossing  angle  between  two  beams  was  10°.  The  crystal  was  oriented  to  have 
its  normal  line  bisecting  the  two  beams.  In  our  experiment  the  transmitted  intensity  of 
the  weak  beam  was  measured.  The  signal  was  detected  by  a  photomultiplier  (  PMT  ) 
and  sent  to  a  602A  DSA  to  be  measured. 

Measurement  and  Data  Presentation 

Experimentally,  we  measured  the  beam  fanning  signal  for  two  different  BSKNN 
crystals  with  no  field  and  with  a  DC  field.  For  the  BSKNN  sample  which  was  heavily 
doped,  the  measured  results  are  shown  in  Fig.4-16  (  no  field  ),  4-17  (  DC  field  3.6 
kv/cm ),  and  4-18  (  DC  field  5.6  kv/cm  ).  The  data  measured  from  Fig.4-16  through 
Fig.4-18  are  listed  in  table  4-1. 
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Table  4-1  Data  of  Beam  Fanning  with  a  DC  Field 


Sample  1 :  BSKNN  heavily  doped 
Wavelength:  458  nm  ( Extraordinary  polarization  ) 


Figure 

Magnitude 

Optical  Transmission 

Response 

No 

of  Field 

Intensity 

Coefficient* 

Time 

Eq  (kv/cm)  I  (mw/cm2) 

T 

x(s) 

Fig.4-16 

0 

250 

3.0% 

3.5 

Fig.4-17 

DC:  3.6 

250 

1.4% 

1.0 

Fig.4-18 

DC:  5.6 

250 

1.1% 

0.46 

*  The  transmission  coefficient  T  is  defined  as  a  ratio  of  the 
transmitted  beam  intensity  at  equilibrium  to  the  initial 
intensity  of  that  beam. 
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Beam  Fanning 
BSKNN  heavily  doped 
I  =  250  mw/cm2 
no  field  ( t  =  3.5  s  ) 


2s/div 

Fig.  4-16 

Note:  This  experimental  curve  was  shown  on  the  screen  of  DSA  oscilloscope.  The 
vertical  axis  gave  the  intensity  of  the  transmitted  beam. 
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Beam  Fanning 
BSKNN  heavily  doped 
I  =  250  mw/cm2 

DC  field:  3.6kv/cm  (xsl.Os) 


- - - - — - — - -  2s/div 

Fig.  4-17 

Beam  Fanning 
BSKNN  heavily  doped 
I  =  250  mw/cm2 

DC  field:  5.6kv/cm  (xs0.46s) 


108mv 

/div 


Fig.  4-18 


2s/div 


Using  an  undoped  BSKNN  sample,  the  measurement  results  are  shown  in  Fig.4-19  (  no 
field  )  and  4-20  (  DC  field  4  kv/cm  ).  For  the  undoped  BSKNN  crystal  we  also 
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measured  a  beam  fanning  signal  while  focusing  the  beam  into  the  crystal  in  order  to 
make  the  optical  intensity  in  the  crystal  about  4  times  stronger.  The  results  are  shown 
in  Fig.4-21  (  no  field )  and  4-22  (  DC  field  4  kv/cm  ).  The  data  measured  from  Fig.4- 
19  through  Fig.4-22  are  listed  in  table  4-2. 

Table  4-2  Data  of  Beam  Fanning  with  DC  Field 
Sample  2:  BSKNN  undoped 


Wavelength:  458  nm  (  Extraordinary  polarization ) 


Figure 

Magnitude 

Optical  Transmission 

Response 

No 

of  Field 

Intensity  Coefficient* 

Time 

Eq  (kv/cm)  I  (mw/cm2)  T 

T(S) 

Fig.4-19  0 

160 

90.0% 

5.6 

Fig.4-20  DC:  4.0 

160 

18.5% 

1.5 

The  following  figures  are  for  beam  fanning  with  focusing. 


Fig.4-21 

0 

640 

90.0% 

1.1 

Fig.4-22 

DC:  4.0 

640 

19.8% 

0.6 

*  The  transmission  coefficient  T  is  defined  as  a  ratio  of  the 
transmitted  beam  intensity  at  equilibrium  to  the  initial 
intensity  of  that  beam. 
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Beam  Fanning 


BSKNN  undoped  1=160  mw/cm2 

no  field  (t  =  5.6s) 


Beam  Fanning  with  Focusing 
BSKNN  undoped 


I  =  640  mw/cm2 
no  field  ( x  =  1.1  s ) 
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As  expected  these  results  show  that  by  focusing  the  speed  of  response  can  be 
decreased.  In  fact,  when  the  spot  size  is  greatly  reduced  to  30  Jim  (  in  SBN  )  the 
response  time  was  extremely  fast.  The  experimental  results  are  shown  in  Fig.4-23, 
which  gave  the  transmitted  intensity  as  a  function  of  time  for  an  input  power  of  about 
0.04  mw.  As  can  be  seen  from  the  curve,  the  characteristic  response  time  is  on  the 
order  of  800  ms. 


Beam  fanning  with  focusing 
Sample:  SBN:60 


Fig.  4-23  Transmitted  Intensity  with  Applying  Field  and  Focusing  Beam 

Comparing  the  data  shown  in  Fig.4-16,  4-17  and  4-18,  we  see  that  with  an 
increase  in  applied  DC  field:  (1)  the  transmission  coefficient  is  decreasing  which 
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means  the  energy  exchange  is  increasing;  and  (2)  the  response  time  is  getting  shorter. 
Similar  results  are  also  obtained  from  the  comparison  of  Fig.4-19  with  Fig.4-20. 

By  comparing  Fig.4-19  with  4-21  and  Fig.4-20  with  4-22,  we  see  that  with 
focusing  the  beam,  the  photorefractive  response  time  gets  shorter  without  a  trade  off  in 
sensitivity.  So  on  the  use  of  an  applied  field,  the  response  time  can  be  further 
decreased  with  focusing  the  laser  beams. 

In  Fig.  4-20  (  also  in  Fig.  4-22  )  we  see  a  field  screening  that  is  caused  by  the 
applied  DC  field.  A  few  seconds  after  turning  the  laser  beam  on,  the  intensity  of  the 
transmitted  beam  starts  increasing.  In  other  words,  the  space-charge  field  was  getting 
smaller.  This  is  due  to  the  electric  field  screening. 

In  Fig.4-23  we  see  an  oscillation  in  the  beam  fanning  signal.  This  agrees  with 
the  theoretical  prediction  in  section  I  of  this  chapter.  According  to  eq.(4-6),  the 
oscillation  frequency  of  the  space-charge  field  is 

©  a  (tr/te)/[  (1+tr/td)2  +  (tr/te)2  ] 

With  no  applied  field  (  1  /xE=0 ),  ©  is  zero.  With  an  increase  in  an  applied  field,  ©  can 
be  seen  to  increase.  In  this  experiment,  we  applied  a  strong  electric  field  (  E0=6  kv/cm 
)  to  the  crystal.  Therefore,  with  a  strong  applied  field  we  can  make  the  frequency  very 
large.  That’s  why  we  see  an  oscillation  in  the  beam  fanning  signal. 

In  the  two  beam  coupling  experiments  we  measured  an  intensity  time  evolution 
of  the  weak  beam  for  the  BSKNN  undoped  crystal.  Fig.4-24  is  the  experimental  curve 
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with  no  field  and  Fig.4-25  is  with  a  DC  field  of  3.6  kv/cm.  The  data  measured  from 
Fig.4-24  and  Fig.4-25  are  listed  in  table  4-3. 


Table  4-3  Data  of  Two  Beam  Coupling  with  DC  Field 


Sample  2:  BSKNN  undoped 
Wavelength:  458  nm  ( Ordinary  polarization  ) 

Full  Crossing  Angle:  20=10° 

Figure  Magnitude  Intensity  Intensity  m  Relative  Response  Modu- 


No  of  Field 

of  Weak 

of  strong 

Energy 

Time 

lation 

Eo 

Beam  Ij 

Beam  I2 

Exchange**  x 

(kv/cm) 

(mw/cm2)  (mw/cm2) 

AIi/I, 

(s) 

Fig.4-24  0 

0.6 

64 

0.2 

0.1 

9.0 

Fig.4-25  DC:  3.6 

0.6 

64 

0.2 

0.9 

5.0 

m=2(I,I2)i/2/(l1+l2) 

*  The  modulation  is  defined  as  a  ratio  of  the  oscillation  amplitude  to 
the  signal  magnitude. 

**  The  relative  energy  exchange  is  defined  as  a  ratio  of  the  intensity 
gain  to  the  original  intensity. 
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The  basic  idea  is  to  examine  if  the  iwo-bcam  coupling  can  be  enhanced  in  both 
magnitude  and  speed  with  an  applied  D  C.  field.  In  fact  wc  will  show  that  for  SBN  and 
BSKNN  this  is  indeed  the  case.  The  problem  however  is  that  screening  of  the  applied 
field  leads  to  a  loss  of  this  enhancement  after  a  time  about  equal  to  the  dielectric  response 
time,  or  about  ten  times  longer  than  the  photorcfiractive  response  time.  Figure  4-24  and  4- 
25  demonstrate  that  the  response  time  and  the  magnitude  of  the  two-beam  coupling  are 
enhanced  with  an  applied  field.  In  fact  the  response  time  improves  from  9  seconds  to  5 
seconds  and  the  energy  transfer  increases  many  fold  over  the  no  field  case.  The  only 
difference  between  figures  4-24  and  4-25  is  that  a  field  of  3.6  KV/cm  is  applied  for  figure 
4-25.  Notice  that  the  enhancement,  however,  is  gradually  lost  as  a  function  of  time  due, 
to  screening  of  the  applied  field. 
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Two  -  beam  coupling  BSKNN  undoped  Total  I  =  64  mw/cm2 

o 

Full  crossing  angle:  10  no  field  (  x  =  9.0  s  ) 


75mv 

/div 


Fig.  4-24  5s/div 


Two  -  beam  coupling  BSKNN  undoped  Total  I  =  64  mw/cm2 

o 

Full  crossing  angle:  10  DC  field:  3.6  kv/cm  ( x  =  5.0  s ) 


We  also  measured  the  signal  of  the  weak  beam  by  focusing  both  beams  into  the  crystal. 
The  results  using  two  focused  beams  are  shown  in  Fig.4-26, 4-27, 4-28  and  4-29.  The 
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data  measured  from  Fig.4-26  through  Fig.4-29  are  listed  in  table  4-4. 

Table  4-4  Data  of  Two  Focused  Beam  Coupling  with  DC  Field 
Sample  2:  BSKNN  undoped 
Wavelength:  458  nm  (  Ordinary  polarization  ) 

Full  Crossing  Angle:  20=10°  (  Focal  Length  of  the  Lens:  100  cm) 

Figure  Magnitude  Intensity  Intensity  m  Relative  Response  Modu- 

No  of  Field  of  Weak  of  strong  Energy  Time  lation* 

Eo  Beam  Beaml2  Exchange**  x 

(kv/cm)  (mw/cm2)  (mw/cm2)  AI1/I1  (s) 


Fig.4-26 

0 

2.6  _ 

260 

0.2 

0.1 

4.0 

Fig.4-27  DC:  3.6 

2.6 

260 

0.2 

1.4 

1.7 

Fig.4-28 

0 

26 

2600 

0.2 

0.3 

0.9 

Fig.4-29  DC:  3.6 

26 

2600 

0.2 

1.3 

0.2 

m=2(I1I2)i'2/(I1+l2) 

*  The  modulation  is  defined  as  a  ratio  of  the  oscillation  amplitude  to 
the  signal  magnitude. 

**  The  relative  energy  exchange  is  defined  as  a  ratio  of  the  intensity 
gain  to  the  original  intensity. 
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Two  focused  beam  coupling  BSKNN  undoped 

o 

Full  crossing  angle:  10  Total  I  =  260  raw/cm2 


ls/div 


DC  field:  3.6  kv/cm 


Fig.  4-26 

(T=  1.7  s) 


BSKNN  undoped 
Total  I  =  2.6  w/cm2 
( x  =  0.9  s ) 


Two  focused  beam  coupling 

o 

Full  crossing  angle:  10 
no  field 


500ms/div 
Fig.  4-28 

DC  field:  3.6kv/cm  (x  =  0.2s) 


79mv 

/div 


500ms/div 
Fig.  4-29 
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By  comparing  Fig.4-24  (  no  field  )  with  Fig.4-25  (  DC  field  ),  Fig.4-26  ( no  field  ) 
with  Fig.4-27  ( DC  field  ),  Fig.4-28  (  no  field  )  with  Fig.4-29  (  DC  field  ),  it  is  shown 
that  an  applied  DC  field  can  really  enhance  the  photorefractive  effect  in  both  the 
magnitude  and  the  speed. 


Direct  Comparison  Of  The  Theoretical  Analysis 
With  The  Experimental  Results 


Experimentally,  we  measured  the  intensity  change  of  the  weak  beam  due  to  two 
beam  coupling.  The  gain  of  the  weak  beam  is  given  by 

lie'll  =  (  1  +  c  )exp(rLcff)/[  1  +  c*exp(rLefr)  ]  ( 4-8 ) 


But  as  the  weak  beam  is  much  weaker  than  the  strong  beam,  then  the  diffraction  signal 
due  to  the  strong  beam  can  not  be  neglected.  Under  this  circumstance,  the  intensity 
measured  from  the  weak  beam  would  include  two  parts:  one  is  the  energy  gain  of  the 
weak  beam  and  the  other  is  the  diffraction  signal  of  the  strong  beam  due  to  the  grating. 
If  the  two-beam  energy-coupling  gain  coefficient  of  the  weak  beam  is  T  and  the 
diffraction  signal  from  the  strong  beam  is  D,  T  in  eq.  (4-8)  would  be  replaced  by  (  F2  + 
D2  )1/2. 


I 

I 


l£ 

I 


( 1+  c  )  e 


Vp  +  D2f 


err 


1  +  c  e 


Vp+D2Lel 


'eff 


(4-9) 
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By  using  eq.  (4-9)  we  can  convert  the  measured  intensity  change  to  (  F2  +  D2  )m.  On 
the  other  hand  we  can  convert  the  space-charge  field  obtained  numerically  to  the  gain 
coefficient  T  of  the  weak  beam.  Then  the  comparison  of  two  gains,  experimental  and 
theoretical,  will  give  us  a  direct  comparison. 

A.  Gain  coefficient  obtained  from  the  experimental  data 

The  experimental  data  we  used  for  deriving  (  F2  +  D2  )1/2are  shown  in  Fig.  4-28 
and  Fig.  4-29.  For  the  data  we  used,  c  is  1/100  and  Lcff  is  about  4  mm.  For  small  gain, 
therefore  according  to  eq.  (4-8),  we  can  get 

Ilc/Ii=exPWp+D2Lerr) 

Then  Vp  +  D2=ln(Ilc/ 1,)  /LelT 

This  equation  enables  us  to  obtain  the  total  magnitude  of  both  the  gain  coefficient 
and  the  diffraction  signal  from  the  intensity  change  which  we  measured 
experimentally.  In  the  experiment,  after  a  grating  was  formed  inside  the  crystal  we 
blocked  the  weak  detected  beam  and  measured  only  the  diffraction  signal  from  the 
strong  beam.  The  diffraction  signal  was  found  to  be  one  fourth  of  the  energy  exchange 
signal.  So  we  can  get  (  F2  +  D2  )1/2  =  ( r2  +  P/4  )  1/2 =  1 . 1 T.  Then  based  on  the  data  in 
Fig.  4-28  (  no  field  )  and  Fig.  4-29  (  DC  field  )  we  calculated  the  corresponding  T 
versus  time,  which  are  shown  in  Fig.  4-30  and  4-32.  From  these  figures  we  can 

determine  the  response  time  for  each  case. 

B.  Gain  coefficient  obtained  from  the  numerical  computation 

Numerically,  by  using  the  same  data  as  in  the  experiment  and  varying  the  parameters 
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Using  a  Cylindrical  Lens  and  an  Applied  Electric  Field. 

While  focusing  coupled  with  an  applied  D.C.  seems  promising  it  suffers  from 
the  screening  effect  seen  in  figures  4-20  and  4-22.  That  is,  alter  some  time  charges 
move  to  the  edges  of  the  incident  beam  and  produce  a  screening  field  which  brings  to 
zero  the  enhancement  of  the  applied  field.  This  can  be  avoid  using  a  c  focus  the 
incident  beam  in  one  direction  only  and  maintains  an  optically  induced  electrical  path 
between  electrodes.  As  a  result,  using  both  an  applied  D.C.  field  along  with  the 
cylindrical  lens  for  focusing  an  enhancement  of  four  orders  of  magnitude  is  possible. 
While  this  is  still  marginal  it  does  accomplish  the  speed  enhancement  without  loss  of 
O.D.  due  to  screening. 

Conclusions 

Our  results,  both  theory  and  experiment,  clearly  show  that  the  application  of  an 
electric  field  can  enhance  both  the  speed  and  the  magnitude  of  the  photorefractive  effect. 
In  fact,  the  data  and  modeling  clearly  show  that  the  use  of  a  sped  cal  lens  to  focus  to 
about  100  micron  spot  size  with  an  applied  field  gives  the  four  order  of  magnitude 
improvement  in  the  response  time  that  we  were  seeking.  This  is  accomplished  without 
sacrificing  O.D.  The  problem  is  that  focusing  with  an  applied  field  leads  to  the 
screening  effect.  However,  the  case  of  an  applied  field  with  a  cylindrical  lens  is  very 
successful  since  the  incident  light  fills  the  crystal. 
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V.  APPLIED  A.C.  ELECTRIC  FIELD 

The  photorefractive  effect  can  be  enhanced  by  applying  a  DC  field  across  a 
crystal.  An  applied  DC  electric  field  can  (1)  enhance  the  magnitude  of  the  space-charge 
field  and  corresponding  photorefractive-index  change  and  (2)  speed  up  the 
photorefractive  response  time.  However,  the  technique  of  applying  a  DC  field  has 
some  drawbacks  or  disadvantages  of  its  own  (Ref.  [l]-[4].  Equation  (4-4)  showed  a 
phase  shift  between  the  space-charge  field  and  the  intensity  pattern.  When  there  is  no 
applied  field,  the  phase  shift  is  90°,  which  gives  the  maximum  energy  exchange 
between  two  beams.  But,  if  a  DC  field  E,,  is  turned  on  the  phase  will  no  longer 
preserve  the  90°  shift  and  energy  exchanged  between  two  beams  is  diminished.  A 
second  problem  is  that  screening  of  the  DC  field  is  observed.  That  is,  when  a  DC  field 
is  applied  across  a  crystal  the  electrons  in  the  conduction  band  move  in  one  direction. 
As  electrons  move  to  the  edge  of  the  illuminated  area  the  separated  positive  and  negative 
charges  will  be  frozen  there  and  a  screening  of  the  applied  field  will  occur.  As  a  result 
of  screening  the  enhancement  in  the  magnitude  and  speed  of  the  photorefractive  effect  is 
lost.  Experimentally,  it  is  easy  to  observe  the  screening  effect.  Fig.  4-20,  Fig.  4-22 
and  Fig.  4-25  in  last  chapter  demonstrates  a  field  screening  with  an  applied  field.  After 
a  short  time  the  enhancement  in  magnitude  of  the  photorefractive  effect  is  lost.  These 
two  disadvantages,  however,  can  be  avoided  using  an  applied  AC  field. 


SCREENING  OF  EXTERNAL  HELD 
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Fig.  4-31  Theoretical  curve 
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Analytical  Solution  For  Space-Charge  Field  With  An 
Applied  Square-Wave  AC  Field  Of  A  Short  Period 


The  applied  square-wave  AC  field  is 

Eac=  E0forO<t<T/2  or 
=  -E0  for  T/2  <  t  <  T 

where  T  is  the  period  of  square-wave  field  and  the  range  of  T  is  xR  «  T  «  x  (  x:  the 
photorefractive  response  time ). 

The  space-charge  field  with  a  DC  field  is  given  by  Eq.  ( 4-2  ) 

E,  ( t )  =  Esc  ( -e icp  •  e  ^  •  e  +  e iip ) 
where  Escis  given  by  Eq.  ( 4-3  ). 

We  define  a  non-dimensional  physical  quantity  E,’  ( t )  as 

E1(t)  =  Esc-E,’(t) 


Then  E,’  at  the  end  of  each  negative  half  period,  i.e.,  when  t  =  nT,  is  given  by  (  see 
Appendix  II ) 


.  -nT/i 


r;  *  ( nT\  1  -e  ><P  -T/t  i<p  io)T/2  -T/2t 

Ei^nT)=( - — )(-e  -e  +  e  -e  -e 


1  -e 


-T/t 


-i<P  -icp  icoT/2  -T/  2t, 

-e  +  e  -e  -e  ) 


(5-1) 


(  n  =  1,  2,  3...  is  the  period  order  of  the  square-wave  field  ) 

Meantime,  E,’  at  the  end  of  each  positive  half  period,  i.e.  when  t  =  (  n  +  1/2  )  T,  is 
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H,’[(n+i)T]=(  '  ~  )•(  e-,%"''-e-"'.e-toT,2.e'T'J') 

2  1  -  e 

-( n+l)T/t 


1  -nT/  t 

1  ~  e  x  ,  -i<P  -T/t 

*  (  e  •  e  -  e  “T-  e 


1  -  e 
1-e 


-T/t 


,  /  1  e  \  /  i9  icp  -  icoT/2  -T/2t  \  f  K  n  \ 

+  ( - ^77- )  (  e  -e  •  e  -e  )  (  j-4) 


1-e 


Since  the  solution  to  the  differential  equation  of  space-charge  field  ( see  eq.  (4-1) )  can 
also  be  written  as 

E,  ( t )  =  Esc  {  [  E/  ( t  o )  -  e i<f> ]  •  e  '( *' l°)/T  •  e  'im(t '  ^  +  e i<fl } 

Therefore,  E,’  during  each  positive  half  period  is 

Ej’  ( t )  =  [  E,’  (  nT)  -  e i<p]  •  e  ',/T  •  e  ’iwl  +  e  i<(>  (0<t’<T/2)  (5-3) 

During  each  negative  half  period 

E,’  ( t )  =  {  E,’  [  n  +  1/2 ) T  ]  +  e"i<p }  •  e’171-  e-ia,t’-  ei<p  ( 0 <  t’  <T/2 )  (  5-4 ) 

Eqs.  (  5-1 ),  ( 5-2 ),  (  5-3  )  and  (  5-4  )  constitute  the  basic  equations  for  space-charge 
field  under  an  applied  square-wave  AC  field.  In  the  case  of  T  «  x  two  results  can  be 
derived. 

1.  In  general  x^»xR,  so  that  tr/td=0.  In  addition,  for  not  very  large  applied 
field  tE»tR,  so  that  xR/xE=0.  We,  therefore,  can  write  the  following 
expression. 
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(5-5) 


2. 


Then  we  get 


Let’s  go  back  to  Ch.4  and  look  at  Eq.  ( 4-5  )  and  Eq.  ( 4-6  ).  It’s  not 
difficult  to  find  that  Eq.  (  5-5  )  yields  a  relation  of  1/x  >  co. 
Therefore,  if  T  «  x,  then  T  «  1/co  or  ooT«l .  This  relation  will  be  used  in 
the  following  derivation. 

In  the  case  of  T  «  x  only  after  many  many  periods  of  AC  field  ( in  other 
words.,  n  must  be  a  very  large  number  )  an  equilibrium  can  be  reached.  So 
we  will  have 


(n+l)T/x  =  nT/x  forn»l 
Thus  Eq.  (  5-2 )  becomes 

1  -  1  _  -nT/x 

E1’[(n  +  ~)  T]  =  (~~~ . 6  -  )•(  e”  e  T/  T-  e~  "*•  e- T/  2-  e‘T/  2  T 

2  1  -  e 

i<P  i<p  -iuT/2  -T/2t  )  (5-6) 

+  e  — e  ■  e  -e  1 


From  either  Eq.  (  5-1  )  or  Eq.  (  5-6  ),  both  will  have  the  same  result  for  the 
imaginary  part 

ImE,  [  nTor  ( n+1/2  )T  ]  = 

(  - — — )  [-eT/Tsin<P  +  e‘T/2‘t  sin  (cot/2  +  cp) 

1  -e’T/T 

+  sin9  +  e'T/2Tsin(coT/2-9)l  (5-7) 

Eq.  (  5-7  )  shows  that  in  the  end  of  each  half  period  of  AC  field  ( whatever 
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positive  or  negative  half  period  )  the  space-charge  field  stands  the  same. 
From  this  analysis  we  can  examine  the  response  time  and  the  phase  shift. 
When  an  equilibrium  state  is  established  n  must  be  a  large  number,  or  in  other 
words,  n  — >  °°  in  Eq.  (  5-7  ),  then  the  equilibrium  value  of  the  imaginary  part 
of  the  space-charge  field  is 

ImEj’  eq  =( - [-eT/tsin<P  +  e'T/2x  sin(G)T/2  +  9) 

1-e 

+  sin9  +  e‘T/2xsin(coT/2-9)]  (5-8) 

When  Im  E,’ eq  reaches  its  e"1  value  or  when  nT/x  =  1,  Eq.  (  5-7  )  becomes 

Im  E,’(nT/x  =  1 )  =  (— — ^— )  [-e‘T/%in9  +  e‘T/2x  sin  (cot/2  +  9) 

1-e  x  ' 

+  sin  9  +  e‘T/2x  sin  (C0T/2  -  9)] 

=  (1-1/e)  ImE,’eq 


(5-9) 

That  is  the  field  reaches  its  e’1  value  after  time  nT/x  =  1  or  nT  =  x.  But  x  is 
also  the  time  for  the  space-charge  field  to  reach  its  e'1  value  for  same  applied 
DC  field.  Therefore,  this  tells  us  that  the  response  time,  in  the  case  of 
applying  a  square-wave  AC  field,  with  a  period  of  T  «  x,  is  the  same  as  in 
the  case  of  applying  a  DC  field. 

To  see  the  phase  shift  in  equilibrium,  for  T  «  x  we  can  make  an 
approximation 


64 


Then  from  Eq.  (5-1  ),  for  n  -+  °° 

*  1 

ReE,  (nT)=( - — )  [_cos  (p+cos((DT/2  +cp) 

1-  e 

-cos  <P+cos(coT/2  +9)1 

1  ,  2 

="( - T77 )  coscp (  ~)  since  coT/2«l 

1-  e  2 

Due  to  the  same  reason,  for  n  ->  «> 

ImE,’  (nT)=  ( — )  Cos(p(  ©T ) 

1-e 

If  we  write  the  complex  number  E,’  (nT)  as 
E,’  (nT)  =1 E,’  (nT)  l-eid 

then  tanf)  =  ImE,’  (nT)/ReE,’  (nT)  =  -4/(coT)  -+  -°°  for  0)T«1 
we  will  get 

■d  =  90°+e  (5-10) 

where  e  is  an  infinitesimal. 

Also  from  eq  (  5-2  ),  by  using  the  same  way  as  above  we  can  get  for  n— 

ReE,’  ( (n+1/2)  T)  =(  — ^ )  coscp(  “!)' 

1-e  2 

ImE,’  ( (n+1/2)  T)  =  ( - ^ )  costp(  (oT  ) 

1-e 

We  write  E,’  ( (n+1/2)  T)  =  IE,’  ( (n+1/2)  T)l  eiv 

tan\[/=  ImE,’  ( (n+1/2)  T)/ReE,’  ( (n+1/2)  T)  =2/(coT)  -+oo 
for  0)T«1 

\|/  =  90°-e  (5-11) 
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Eqs  (  5-10  )  &  (  5-1 1  )  show  that  the  phase  shift  is  always  about  90°  in  both  the 
positive  and  negative  half  periods.  Physically,  this  means  that  we  can  get  a 
maximum  energy  exchange  between  two  beams  by  applying  a  square- wave  short 
period  AC  field  to  the  crystal.  This  is  a  remarkable  result. 

Numerical  Solution  For  Space-Charge  Field  With  An 
Applied  Square  Wave  AC  Field  Of  An  Arbitrary  Period 

We  have  already  obtained  an  analytical  solution  for  the  space-charge  field  with  an 
applied  square  wave  AC  field  of  a  very  short  period  (  T«t  ).  But,  in  general,  the 
photorefractive  effect  can  be  enhanced  by  applying  an  AC  field  with  any  period. 
Therefore,  we  have  to  resort  to  numerical  solution  for  the  space-charge  field  with  an 
external  square  wave  AC  field  of  an  arbitrary  period.  Considering  that  the  applied  AC 
field  changes  its  value  through  a  process  of  +E0— >0— >-Eq  therefore,  we  must  make  a 
transition  between  the  positive  and  the  negative  fields  when  we  look  for  a  numerical 
solution  for  the  space-charge  field. 

The  crystal  physical  parameters  for  the  AC  field  case  are  same  as  those  for  the 
DC  field  case  in  Chapter  4.  The  only  difference  is  that  a  fixed  DC  field  value  Eo  is 

now  replaced  by  an  alternating  field  value  +Eo.  Fig.5-1,  5-2  and  5-3  show  the  space- 
charge  field  time  evolution  with  an  AC  field  for  different  NA.  Fig.5-4  and  5-5  show 
the  space-charge  field  with  an  AC  field  of  different  frequency. 
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The  numerical  solutions  for  an  AC  field  tell  us  the  following  points. 

1.  General  effect  of  AC  field 

An  applied  AC  field  can  give  an  efficient  enhancement  on  the  photorefractive 
effect  in  both  magnitude  and  speed.  This  can  be  seen  by  comparing  Fig.5-1  (  AC 
field  )  with  Fig.4-1  (  no  field  ),or  Fig.5-2  (  AC  field  )  with  Fig.4-3  (  no  field  ),or 
Fig.5-3  (  AC  field  )  with  Fig.4-5  ( no  field  ).  Also,  by  comparing  Fig.5-1  (  AC  field 
)  with  Fig.4-2  (  DC  field  ),  or  Fig.5-2  (  AC  field  )  with  Fig.4-4  (  DC  field  ),  or 
Fig.5-3  (  AC  field  )  with  Fig.4-6  (  DC  field  ),  we  can  see  that  the  effect  of  AC  field 
on  enhancing  the  magnitude  and  response  time  of  the  space-charge  field  is  not  as 
strong  as  that  of  DC  field,  which  can  be  explained  physically.  From  the  analytical 
solution  presented  in  last  section,  the  equilibrium  value  of  the  space-charge  field 
with  an  AC  field  is  given  by  eq.  (5-8),  from  which  we  can  get 
ImE, cq  =  sincp  +  2-[l /(eT/2x  -e  't/2t  )]-sin(coT/2)-cos(p 
In  this  equation,  when  T  increases  the  equilibrium  value  of  the  space-charge  field 
decreases.  As  stated  in  last  section,  when  T«t  the  space-charge  field  with  an  AC 
field  is  almost  same  as  that  with  a  DC  field.  So  when  T  is  not  very  small  compared 
to  t,  the  space-charge  field  will  be  less  than  the  case  of  T«x  This  means  that  the 
space-charge  field  with  an  AC  field  having  a  not  very  short  period  is  not  as  strong 
as  that  with  a  DC  field.  That’s  why  the  effect  of  AC  field  on  enhancing  the 
magnitude  of  the  space-charge  field  is  not  as  strong  as  that  of  DC  field.  As  to  the 
response  time,  the  physical  explanations  are  as  follows.  During  the  first  half  period 
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of  an  AC  field,  the  temporal  behavior  of  the  space-charge  field  is  same  as  the  DC 
curve.  But  immediately  after  the  first  half  period,  the  space-charge  field  intends  to 
decrease  and  an  oscillation  starts.  The  oscillation  makes  the  response  time  with  AC 
field  not  as  fast  as  that  with  DC  field. 

2.  Influence  of  the  trap  number  density 

The  numerical  solutions  for  different  trap  number  density  of  the  crystal  ( Fig.5-1,  5- 
2  and  5-3  )  tell  us  that  in  order  to  get  a  stronger  enhancement  of  the  photorefractive 
effect  with  an  applied  AC  field  it  is  important  that  a  crystal  must  be  doped  to  have 
appropriate  donor  number  density  and  trap  number  density  as  in  the  DC  case.  In 
these  three  figures,  the  peak  performance  for  the  response  time  is  that  with  NA 
value  of  1.5*1016  cm-3.  The  physical  reason  is  same  as  what  we  discussed 
previously  in  DC  field  case  in  Chapter  4. 

3.  Effect  of  AC  field  frequency 

All  these  figures  ( Fig.5-1  through  Fig.5-5  )  show  that  the  frequency  of  total  space- 
charger  field  is  twice  the  frequency  of  an  applied  AC  field.  The  numerical  results 
with  different  frequency  of  an  AC  field  (  Fig.5-4  and  5-5  )  indicate  that  an  increase 
of  an  AC  field  frequency  will  cause  less  modulation  in  the  space-charge  field  and  a 
faster  response  time.  As  stated  previously,  with  an  AC  field  the  space-charge  field 
would  intend  to  decrease  and  start  to  oscillate  immediately  after  the  first  half  period 
of  the  AC  field.  With  a  low  frequency  of  AC  field  the  oscillation  gets  stronger  than 
that  with  a  high  frequency.  That’s  why  a  high  AC  field  frequency  causes  less 


modulation  in  the  space-charge  field  and  a  faster  response  time. 

Experimental  Setup  And  Data  Analysis  For  AC  Field 

The  experimental  apparatus  using  an  AC  field  is  mostly  same  as  that  using  a  DC 
field.  The  AC  voltage  was  produced  by  using  HEWLETT-PACKARD  8002A  PULSE 
GENERATOR,  which  was  then  amplified  1000  times  using  a  610C  H.V.  SUPPLY 
AMPLIFIER.  The  amplified  AC  high  voltage  was  then  connected  to  the  crystal.  We 
investigated  the  effect  of  AC  field  on  the  photorefractive  phenomena  through  both 
beam  fanning  and  two  beam  coupling  experiments.  The  experimental  setup  is  same  as 
in  Fig.4-14  (  for  beam  fanning )  and  Fig.4-15  (  for  two  beam  coupling ). 

Experimental  Data 

The  beam  fanning  curves  are  shown  in  Fig.5-6  through  Fig.5-1 1  for  a  BSKNN 
heavily  doped  crystal,  Fig.5-12  through  Fig.5-14  for  a  BSKNN  undoped  crystal,  and 
Fig.5-15  and  Fig.5-16  for  a  BSKNN  undoped  crystal  using  a  focused  laser  beam.  The 
data  measured  from  these  figures  are  listed  in  table  5-1. 
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Table  5- 1  Data  of  Beam  Fanning  with  AC  Field 


Sample  1:  BSKNN  heavily  doped 
Wavelength:  458  nm  ( Extraordinary  polarization  ) 

Figure  Magnitude  Frequency  Optical  Transmission  Response 


No  of  Field 

Eq  (kv/cm) 

of  Field 

f  (Hz) 

Intensity  Coefficient* 

I  (mw/cm2)  T 

Time 

x(s) 

Fig.5-6  0 

0 

250  3.0% 

3.5 

Fig.5-7  AC:  3.6 

0.33 

250  6.4%(top)  3.0%(bottom)  2.7 

Fig.5-8  AC:  3.6 

25 

250  2.6%(top)  1.5%(bottom)  2.6 

Fig.5-9  AC:  5.6 

0.33 

250  5.5%(top)  2.5%(bottom)  2.3 

Fig.5-10  AC:  5.6 

25 

250  2.2%(top)  1.4%(bottom)  2.3 

Fig.5-1 1  AC:  5.6 

50 

250  2.1%(top)  0.8%(bottom)  2.3 

*  The  transmission  coefficient  T  is  defined  as  a  ratio  of  the 
transmitted  beam  intensity  at  equilibrium  to  the  initial 
intensity  of  that  beam. 


72 


Beam  Fanning  BSKNN  heavily  doped  I  =  250  mw/cm2 
AC  field:  5.6  kv/cm  0.33  Hz  ( x  =  2.3  s ) 


lOOmv 

/div 


Fig.  5-9 

AC  field:  5.6  kv/cm  25  Hz  ( x  =  2.3  s  ) 

( 


81mv 


2s/div 
Fig.  5-10 

AC  field:  5.6  kv/cm  50  Hz  ( x  =  2.3  s  ) 
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Fig.  5-11 
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By  comparing  Fig.5-6  (  no  field  )  with  Fig.5-7  through  Fig.5-11  (  AC  field  )  we 
see  that  the  photorefractive  effect  in  both  magnitude  and  speed  is  enhanced  by  applying 
an  AC  field.  But  from  comparing  Fig.5-7  and  Fig.5-8  (  AC  field  )  with  Fig.4-17  (  DC 
field ),  or  Fig.5-9  through  Fig.5-1 1  (  AC  field  )  with  Fig.4-18  ( DC  field  ),  it  shows  that 
the  effect  of  AC  field  on  enhancing  the  magnitude  and  speed  of  the  photorefractive 
effect  is  not  as  strong  as  that  of  DC  Field.  A  detailed  comparison  of  Fig.5-7  (  AC  field 
of  3.6  kv/cm )  with  Fig.4-17  ( DC  field  of  3.6  kv/cm ),  Fig.5-9  (  AC  field  of  5.6  kv/cm 
)  with  Fig.4-18  (  DC  field  of  5.6  kv/cm  )  showed  that  the  fanning  signal  in  the  first  half 
period  of  AC  field  was  almost  same  as  DC  figure,  but  after  that  an  oscillation  happened 
in  the  fanning  signal  with  AC  field.  The  oscillation  then  made  the  response  time  longer 
than  that  of  DC  field  case. 

By  comparing  Fig.5-7  and  Fig.5-8  (  AC  field  of  3.6  kv/cm  )  with  Fig.5-9  through 
Fig.5-1 1  ( AC  field  of  5.6  kv/cm  ),  it  can  be  seen  that  an  increase  of  AC  field  amplitude 
will  cause  a  stronger  enhancement  of  the  photorefractive  effect  in  both  magnitude  and 
speed.  Fig.5-9  (  or  Fig.5-10  )  with  AC  field  of  5.6  kv/cm  showed  a  stronger  beam 
fanning  signal  and  less  response  time  than  Fig.5-7  (  or  Fig.5-8  )  with  AC  field  of  3.6 
kv/cm. 

By  comparing  Fig.5-7  (  AC  field  of  0.33  Hz )  with  Fig.5-8  (  AC  field  of  25  Hz  ), 
or  Fig.5-9  (  AC  field  of  0.33  Hz  ),  Fig.5-10  (  AC  field  of  25  Hz  )  and  Fig.5-11  (  AC 
field  of  50  Hz )  with  each  other,  it  can  be  found  that  an  increase  of  AC  field  frequency 
will  reduce  the  modulation  of  oscillation  in  signal. 

All  these  experimental  results  are  in  general  agreement  with  the  numerical 
solutions  presented  in  last  section  (  Section  II  of  this  chapter ).  Another  beam  fanning 
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experiment  was  performed  for  the  undoped  BSKNN  sample.  The  measurement  results 
are  shown  in  Fig.5-12  through  Fig.5-14  for  beam  unfocused,  and  Fig.5-15  and  Fig.5-16 
for  beam  focused.  The  data  measured  from  these  figures  are  listed  in  table  5-2. 

Table  5-2  Data  of  Beam  Fanning  with  AC  Field 
Sample  2:  BSKNN  undoped 
Wavelength:  458  nm  ( Extraordinary  polarization  ) 

Figure  Magnitude  Frequency  Optical  Transmission  Response 
No  of  Field  of  Field  Intensity  Coefficient*  Time 

Eo  (kv/cm)  f  (Hz)  I  (mw/cm2)  T  x  (s) 


Fig.5-12  0  0  160  90.0%  5.6 

Fig.5-13  AC:  4.0  _  0.4  160  39.1%(top)  17.5%(bottom)  4.0 
Fig.5-14  AC:  4.0  20  160  83.2%(top)  79.8%(bottom)  3.8 


The  following  figures  are  for  beam  fanning  with  focusing. 


Fig.5-15  0  0  640  90.0%  1.1 

Fig.5-16  AC:  4.0  1.2  640  60.7%(top)  41.8%(bottom)1.0 


*  The  transmission  coefficient  T  is  defined  as  a  ratio  of  the  transmitted 
beam  intensity  at  equilibrium  to  the  initial  intensity  of  that  beam. 
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Beam  Fanning  BSKNN  undoped  I  =  160  mw/cm2 
no  field  ( x  =  5.6  s ) 


2s/div 
Fig.  5-14 
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By  comparing  Fig.5-12  (  no  field  )  with  Fig.5-13  and  Fig.5-14  (  AC  field  ),  or 
Fig.5-15  (  no  field )  with  Fig.5-16  (  AC  field  ),  it  can  be  seen  that  an  applied  AC  field 
really  enhance  the  photorefractive  effect  in  both  magnitude  and  speed. 

By  comparing  Fig.5-13  and  Fig.5-14  (  AC  field  )  with  Fig.4-20  (  DC  field  ),  or 
Fig.5-16  ( AC  field  and  focusing  )  with  Fig.4-22  ( DC  field  and  focusing  ),  we  see  that 
although  the  effect  of  AC  field  on  enhancing  the  photorefractive  effect  is  not  as  strong 
as  that  of  DC  field,  but  with  an  applied  AC  field  the  field  screening  is  avoided,  which 
happens  in  the  case  of  DC  field  shown  in  Fig.4-20  and  Fig.4-22.  Let’s  compare  Fig.5- 
16  ( AC  field  of  4  kv/cm  )  with  Fig.4-22  ( DC  field  of  4  kv/cm  ).  In  Fig.5-16  the  beam 
fanning  signal  during  the  first  half  period  of  AC  field  is  almost  the  same  as  DC  figure 
shown  in  Fig.4-22.  Then,  in  Fig.4-22  the  signal  keeps  going  down,  but  in  Fig.5-16 
after  the  first  half  period  the  signal  goes  up  and  an  oscillation  happens.  Due  to  the 
oscillation  the  response  time  gets  longer  in  AC  field  case  than  that  in  DC  field  case. 
That’s  why  the  effect  of  AC  field  on  enhancing  the  photorefractive  effect  is  not  as 
strong  as  that  of  DC  field. 

In  Fig.5-13  we  see  that  the  frequency  of  the  beam  fanning  signal  is  twice  the 
frequency  of  an  applied  AC  field,  which  is  in  agreement  with  the  theoretical  prediction 
in  last  section. 

From  comparing  Fig.5-13  and  Fig.5-14  ( beam  unfocused  )  with  Fig.5-16  (  beam 
focused  ),  it  can  be  seen  that  on  the  use  of  AC  field,  if  laser  beam  is  focused,  the 
response  time  can  be  further  decreased. 
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For  the  BSKNN  undoped  crystal,  the  experimental  results  of  two  beam  coupling 
are  shown  in  Fig.5-17  and  Fig.5-18  for  the  case  of  beam  unfocused,  and  Fig.5-19 
through  Fig.  5-24  for  the  case  of  beam  focused.  The  data  measured  from  Fig.5-17  and 
Fig.5-18  are  listed  in  table  5-3,  and  data  from  Fig.5-19  through  Fig.5-24  are  in  table  5- 
4. 

Table  5-3  Data  of  Two  Beam  Coupling  with  AC  Field 
Sample  2:  BSKNN  undoped 
Wavelength:  458  nm  ( Extraordinary  polarization  ) 

Full  Crossing  Angle:  20=10° 

Figure  Magnitude  Frequency  Intensity  Intensity  m  Relative  Response  Modu- 
No  of  Field  of  Field  of  Weak  of  strong  Energy  Time  lation* 

Eq  f  Beam  Ij  Beam  I2  Exchange**  x 

(kv/cm)  (Hz)  (mw/cm2)  (mw/cm2)  AIj/I,  (s) 


Fig.5-17  0  0 

0.6 

64 

0.2 

0.1 

9.0 

Fig.5-18  AC:  3.6  1.2 

0.6 

64 

0.2 

0.9 

6.5 

m  =2(  IjL*)!'2/  (I1+I2) 

*  The  modulation  is  defined  as  a  ratio  of  the  oscillation  amplitude  to  the  signal 
magnitude.  **  The  relative  energy  exchange  is  defined  as  a  ratio  of  the 
intensity  gain  to  the  original  intensity. 
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Two  -  beam  coupling  BSKNN  undoped  Full  crossing  angle:  10 
no  field  Total  I  =  64  mw/cm2  ( x  =  9.0  s ) 


5s/div 
Fig.  5-17 


Two  -  beam  coupling  BSKNN  undoped  Full  crossing  angle:  10 
AC  field:  3.6kv/cm  1.2  Hz  Total  I  =  64  mw/cm2  (is  6.5  s) 


5s/div 
Fig.  5-18 
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Table  5-4  Data  of  Two  Focused  Beam  Coupling  with  AC  Field 
Sample  2:  BSKNN  undoped 
Wavelength:  458  nm  ( Extraordinary  polarization  ) 

Full  Crossing  Angle:  20=10° 

( Focal  Length  of  the  Lens:  100  cm) 

Figure  Magnitude  Frequency  Intensity  Intensity  m  Relative  Response  Modu- 
No  of  Field  of  Field  of  Weak  of  strong  Energy  Time  lation* 
Eq  f  Beam  I,  Beam  I2  Exchange**  x 

(kv/cm)  (Hz)  (mw/cm2)  (mw/cm2)  AI,/I,  (s) 


Fig.5-19  0 

0 

2.6 

260 

0.2 

0.1 

4.0 

Fig.5-20  AC:  3.6 

1.2 

2.6 

260 

0.2 

1.3 

1.7 

0.42 

Fig.5-21  AC:  3.6 

1.2 

2.6 

410 

0.16 

0.9 

1.26 

1.1 

Fig.5-22  AC:  3.6 

1.8 

2.6 

410 

0.16 

1.2 

0.88 

0.55 

Fig.5-23  0 

0 

26 

2600 

0.2 

0.3 

0.9 

Fig.5-24  AC:  3.6 

3.2 

26 

2600 

0.2 

1.7 

0.3 

0.7 

m=2(  1,12)1/2/ (I1+I2) 

*  The  modulation  is  defined  as  a  ratio  of  the  oscillation  amplitude  to  the 
signal  magnitude. 


**  The  relative  energy  exchange  is  defined  as  a  ratio  of  the  intensity  gain  to 
the  original  intensity. 
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By  comparing  Fig.5-18  (  AC  field  )  with  Fig.5-17  (  no  field  ),  or  Fig.5-20  (  AC 
field  )  with  Fig.5-19  ( no  field  ),  or  Fig.5-24  (  AC  field  )  with  Fig.5-23  (  no  field  ),  we 
see  that  an  AC  field  enhances  the  photorefractive  effect  in  both  the  magnitude  and 
speed.  Fig.5-17,  Fig.5-19  and  Fig.5-23  show  almost  no  or  very  little  energy  exchange 
between  two  beams.  But  with  applying  an  AC  field,  Fig.5-18,  Fig.5-20  and  Fig.5-24 
show  an  obvious  energy  exchange  between  two  beams  and  also  a  decrease  of  the 
response  time  compared  with  the  no  field  case. 

By  comparing  Fig.5-18  ( AC  field,  x=6.5  s  )  with  Fig.4-25  (  DC  field,  x=5.0  s  ),  or 
Fig.5-20  (  AC  field,  x=1.7  s  )  with  Fig.4-27  (  DC  field,  x=1.7  s  ),  or  Fig.5-24  (  AC  field 
,  x=0.3  s  )  with  Fig.4-29  (  DC  field,  x=0.2  s  ),  we  see  that  the  effect  of  AC  field  on 
enhancing  the  photorefractive  response  time  is  not  as  strong  as  that  of  DC  field,  but  two 
results  are  very  close.  Also,  the  advantage  of  using  AC  field  is  that  with  the  use  of  AC 
field  a  field  screening  is  avoided,  which  can  be  seen  from  comparing  Fig.5-18  (  AC 
field  )  with  Fig.4-25  (  DC  field  ). 

By  comparing  Fig.5-21  ( AC  field  of  1.2  Hz )  with  Fig.5-22  (  AC  field  of  1.8  Hz  ), 
we  see  that  an  increase  of  AC  field  frequency  can  reduce  the  response  time  and 
enhance  the  energy  exchange  between  two  beams. 

By  comparing  Fig.5-20,  Fig.5-21,  Fig.5-22  and  Fig.5-24  (  AC  field  and  beam 
focused  )  with  Fig.5-18  (  AC  field  and  beam  unfocused  ),  we  find  that  if  the  laser 
beams  are  focused  the  response  time  can  be  further  reduced.  The  results  show  that 
focusing  increases  the  speed  of  response  and  that  the  applied  AC  field  makes  possible 
strong  coupling  simultaneously. 
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Direct  Comparison  Of  Experiment  With  Theory 

By  using  the  same  way  and  same  parameters  as  in  Section  IV  of  Chapter  4,  we  can 

obtain  the  gain  coefficient  for  AC  field  case  from  both  the  experimental  and  numerical 
data.  The  results  are  shown  in  Fig.5-25  and  5-26,  which  are  based  on  the  experimental 
data  of  Fig.5-24.  Two  curves  give  almost  the  same  energy  gain  and  same  response 
time. 


Two  focused  beam  coupling 
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Conclusions 


Both  theory  and  experiment  demonstrate  that  the  application  of  an  A.C.  electric 
field  can  be  used  to  enhance  both  the  speed  and  magnitude  of  the  photorefractive  effect 
without  the  draw  back  of  the  screening  phenomena  observed  using  a  D.C.  electric  field. 
Our  experiments  were  carried  out  using  BSIKNN  as  opposed  to  SBN  or  BaTi03  due  to 
the  higher  depoling  field  of  BSKNN.  Our  initial  experiments  using  SBN  were  not 
successful  after  about  lkv/cm.  However,  for  BSKNN  it  was  possible  to  apply  fields 
up  to  6kv/cm  without  significant  depoling.  It  would  be  interesting  to  find  additional 
materials  with  higher  depoling  fields  since  they  would  allow  even  higher  applied  A.C. 
electric  fields  and  furancement. 

The  clear  result  in  this  investigation  is  that  a  focused  laser  beam  into  a  BSKNN 
sample  enhanced  the  response  time.  Our  experiments  show  that  for  an  incident  beam 
5cm  x  5cm  in  size  the  use  of  a  spherical  focusing  lens  to  focus  to  a  50  micron  beam 
diameter  would  increase  the  intensity  by  six  orders  of  magnitude  in  the  crystal  and 
reduce  the  response  factor.  Since  this  is  much  greater  than  the  four  orders  of 
magnitude  required  this  technique  can  act  as  an  effective  limiter. 
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VI.  DETERMINATION  OF  TRAP  NUMBER  DENSITY 

One  of  the  by-products  of  our  investigation  using  applied  D.C.  electric  fields 
has  been  of  a  new  technique  to  determine  the  trap  number  density  in  photorefractive 
materials.  Our  technique  is  based  on  the  energy  exchange  between  two  laser  beams  in  a 
two-beam  coupling  experiment. 

Two  laser  beams  intersect  inside  a  crystal  and  write  an  interference  pattern. 
Then  due  to  the  charge  migration,  a  change  of  the  refractive  index  is  caused  and  an 
index  grating  is  produced  in  the  crystal,  which  is  shifted  in  phase  with  the  intensity 
pattern.  The  phase  shift  between  the  grating  and  the  fringe  pattern,  thus,  causes  an 
energy  exchange  between  two  beams. 

In  a  two  beam  coupling  experiment  we  cam  measure  the  energy  exchange  signal 
between  two  beams,  Ig  and  the  diffraction  signal  due  to  the  grating,  Id  (Ref.  [1]).  The 
energy  exchange  signal  is  the  intensity  change  of  the  detected  beam  after  two  beams  are 
turned  on,  and  the  diffraction  signal  is  the  intensity  measured  immediately  after  the 
detected  beam  blocked.  If  two  beams  have  equal  intensity  of  I0  and  the  crossing  angle 
between  two  beams  are  small,  Ie  and  Id  must  be  much  less  than  I0.  Then  the  phase  shift 
<|)  between  the  grating  and  the  fringe  pattern  can  be  determined  by  the  following 
equation  .  (Se  Appendix  III). 

sin  <(>  =  Ie/[2(I0Id)1/2]  (6-1) 

Energy  Transfer  Measurement  And  Determination  Of  The 

Phase  Shift  Betwen  The  Grating  And  The  Fringe  Pattern 

The  experimental  set  up  for  two  beam  coupling  is  shown  in  Fig.  6-1. 
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splitter 


Fig.  6-1  Experimental  setup  for  measuring  phase-shift 


The  He-Ne  laser  (  5  mw  of  power  and  X  =  632.8  nm )  was  ordinary  polarized  in  order 
to  minimize  beam  fanning  and  resulting  fluctuations  in  the  two  beam  coupling 
signals.  After  the  polarizer  the  laser  beam  was  split  into  two  beams  IT  and  ID  that 
intersect  inside  a  thin  sample  of  SBN  60  with  0.015%  Cesium  as  a  dopant  (  static 
dielectric  constant  e’  =  950 )  such  that  the  grating  vector  is  parallel  to  the  c-axis.  The 
beam  splitter  is  50%  to  50%  to  make  the  beam  IT  having  the  same  intensity  as  ID.  The 
angle  20  between  the  two  beams  was  22°.  The  beam  Lj(  d  )  was  detected  with  a 
photodiode.  For  all  our  observed  intensity  values  the  photodiode  was  linear.  To 
avoid  the  observation  of  multi-exponential  growth  and  decay  times,  the  detected 
beam  was  expanded  with  a  beam  expander  so  that  only  the  uniform  portion  of  the 
beam  was  detected  by  the  photodiode.  The  output  of  the  diode  was  fed  into  the  602A 
Digitizing  Signal  Analyzer.  After  turning  on  both  beams  a  growth  curve  was 
recorded  and  the  energy  exchange  signal  Ie  was  measured.  When  the  grating  reached 
its  steady  state  we  blocked  the  detected  beam.  The  transmitted  beam  immediately 
went  to  zero,  and  the  diffraction  signal  Id  was  detected.  Finally,  the  intensity  of  the 
transmitted  beam  I0  was  measured.  By  substituting  I0,  I£  and  Id  into  Equation  (6-1  ), 
we  can  calculate  the  phase  shift  between  the  grating  and  the  interference  pattern  in  the 
crystal. 

III.  Determination  Of  Trap  Number  Density  By  Phase  Shift  Measurement  And 
Field  Dependence  Of  Dielectric  Constant  And  Electro-Optic  Coefficient 
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As  discussed  in  chapter  4  and  5,  a  trap  number  density  NA  of  a  crystal  plays  a  key 
role  in  the  photorefractive  response  time.  Therefore,  measuring  NA  accurately 
becomes  an  important  topic  in  the  research  on  photorefractive  phenomena.  In  this 
section  we  will  present  a  new  technique  for  the  measurement  of  NA  by  applying  a  DC 
field  across  a  crystal  and  determining  the  resulted  phase  shift  between  the  grating  and 
the  fringe  pattern. 

After  applying  a  DC  field  across  a  crystal  we  carried  out  the  measurement 
described  in  the  last  section  yielding  values  for  I0,  I£  and  Id.  By  substituting  them  into 
Eq.  (6-1  )  we  calculated  both  the  sin  <p  and  the  tan  cp.  In  order  to  get  accurate  results 
we  varied  the  magnitude  of  DC  field  and  measured  I0,  I£  and  Id  for  each  DC  field 
value  Eg.  The  data  for  different  E0  and  the  calculated  sin  cp,  tan  cp  are  listed  in  Table 
6-1. 

From  the  band  transport  theory  we  have 

tan<p=t+(f2+vi  (6-2) 

^=(tan<‘,-^,/(^2+E»)  (6-3) 

Expression  ( 6-2 )  gives  a  value  of  Eg  for  minimum  cp  of  E0=[EDEq+ED2] 1/2 .  Using 
eq.  (  6-3  )  we  are  able  to  calculate  a  set  of  data  of  1/E?  vs  E0  Eq  is  defined  as 


Eq  =  4neNefr/s  K„ 
N. 


where  (  N  rr= - ^ - ) 

eff  1  +  Na/  ( Njj  Na) 
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Table  6-1  Relation  Between  the  Applied  DC 
Field  and  the  Resulted  Phase  Shift 


Eq  ( v/cm ) 

tantp 

sirup 

0.0 

25.4680 

0.9992 

277.78 

5.48783 

0.9838 

555.56 

3.21020 

0.9548 

833.33 

2.66599 

0.9363 

1111.11 

2.46498 

0.9267 

1388.89 

2.26533 

0.9148 

1666.67 

2.41280 

0.9238 

1944.44 

2.53525 

0.9303 

2222.22 

2.71445 

0.9386 

2511.11 

2.72540 

0.9388 

If  Eq  depends  on  E0,  so  the  dielectric  constant  e  must  depend  on  E0.  For  the 
parameters  used  in  our  experiment  we  found 

Kg=  ^L(2sin0)  =  3.79xlO4cm‘1 

KBTKg  y/ 

Ed=  q  ~ =  963.24  /cm 

E„ r^f^(“T!^)=(4-776x  I0’U)('^£L)  (6'4) 

q  Kg  e  E 
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We  assume  that  the  dependence  of  the  dielectric  constant  e  on  the  magnitude  of  DC 
field  Eq  has  a  form  of 

e(Eo)=  e’d  +  aEo+PEo2) 

=  e’[l+f(Eo)]  (6-5) 

where  a  and  (3  are  constants  to  be  determined. 

By  combining  Eqs  ( 6-3  ),  ( 6-4 )  and  ( 6-5  )  we  can  calculate  1  /  Eq  and  Neff  /  e  ’  for 
each  value  of  DC  field  E0.  The  results  are  listed  in  Table  6-2. 


Table  6-2 


i 

(E0)i 

(v/cm) 

d/Eq)i 

(xlO'^cm/v) 

(l/EqV(l/Eq),= 

{ 1  +f[(E0>j] }/( 1  +f  [(E0) !  ] } 

IWe* 

(x  10  ^cm" 

1 

277.78 

5.5838 

1.0000 

4.71 

2 

555.56 

6.6335 

1.1880 

4.53 

3 

833.33 

7.7571 

1.3892 

4.30 

4 

1111.11 

8.2114 

1.4706 

4.37 

5 

1388.89 

7.6419 

1.3685 

4.92 

6 

1666.67 

8.2526 

1.4779 

4.66 

7 

1944.44 

8.4236 

1.5086 

4.57 

8 

2222.22 

8.6410 

1.5475 

4.37 

9 

2511.11 

8.1296 

1.4559 

4.44 

94 


Based  on  the  data  in  Table  6-2  we  plotted  {  1  +  f  [  (  Eq  )  s  ]  }  j  {  1  +  f  [  (  E0 )  j  ]  } 
versus  ( E0)  j  which  is  shown  in  Fig  6-2.  By  fitting  the  data  we  found  that 
a  =  0.0009515  and  B  = -2.51  xlO'7. 


Applied  Field  ( v/cm ) 

Fig.  6-2 

As  a  function  of  an  applied  DC  field,  the  dielectric  constant  €  has  a  form  of 
e(  Eq)  =  €’  [  1  +0.0009515  E0 -2.51  x  10  *7  Eo2] 


llirail 


4.776  -10‘ 


1  +  0.0009515  E0  -  ( 2.51  10  '7)  E02 


For  each  pair  of  E0  and  Eq  we  calculated  the  value  of  Neff/ e which  is  listed  in  the 


middle  column  of  Table  6-3. 
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Table  6-3  Full  crossing  angle:  22° 


Applied  DC  field  E0  ( v/cm )  Neff/<=  ’  ( x  10 13cm'3)  Neff  ( x  10  16cm'3 ) 


277.78 

4.71 

4.47 

555.56 

4.53 

4.30 

833.33 

4.30 

4.09 

1111.11 

4.37 

4.15 

1388.89 

4.92 

4.67 

1666.67 

4.66 

4.43 

1944.44 

4.57 

4.34 

2222.22 

4.37 

4.15 

2511.11 

4.44 

4.22 

*  The  above  listed  results  are  based  on  the  dielectric  constant  e  being  e  = 

€  ’  ( 1  +  9.2  •  10  ^Eq- 2.51  •  IO^Eq2),  where  s’  is  950  for  our  crystal. 

By  taking  the  average,  we  obtained 

Neff/e  ’  =  ( 4.54  ±  0.38  )  x  10 13  ( cm'3 ). 

For  our  crystal  e  ’  is  950.  Therefore,  we  get  Neff  =  ( 4.3 1  ±  0.36 )  x  10  16  ( cm'3 ). 

Neff  is  approximately  the  trap  density  NA  for  our  crystal.  So  we  can  get 

Na  =  ( 4.31  +  0.36 )  x  10  16  ( cm'3 ).  In  Ref.[3]  NA  value  of  ( 4.2  +  0.2 )  x  10  16  (  cm'3 

)  was  reported  for  a  crystal  of  SBN:60  with  0.015%  Rh.  concentration.  Our  results 
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are  very  closer  to  theirs.  We  plotted  curves  of  sin  cp  vs  E0,  tan  (p  vs  E0  and  e  vs  E, 
respectively.  They  are  shown  in  Fig  6-3,  6-4  and  6-5. 


Beside  the  information  of  the  electric  field  dependence  of  the  dielectric 
constant  which  we  obtained  from  this  experiment,  from  the  measurement  of 
diffraction  signal  we  are  also  able  to  determine  the  electric  field  dependence  of  the 
electro-optic  coefficient.  The  diffraction  efficiency  t|  is  given  by  ( Ref.  [1],  [2] ) 

Tj  =  sin2(  Kd ) 


where  k  =  n  •  An  /  X  cosO  and  An  =  ( 1/2 )  n03  reff  Esc-  m. 

From  the  above  equation  we  can  get  an  expression  for  electro-optic  coefficient  reff 
which  is 


reff  -(• 


sirfVr K  , 
E  n 


2Xcos0, 

ttn^dm 


As  we  applied  a  DC  field  Eq  across  a  crystal,  the  diffraction  signal  T|  and  the  space- 
charge  field  Esc  will  change  with  Eq  but  keep  all  other  parameters  unchanged. 
Therefore,  by  calculating  T|  and  Esc  we  are  able  to  determine  the  dependance  of 
electro-optic  coefficient  reff  on  the  applied  DC  field  E0.  Since  the  laser  beams  are  in 
ordinary  polarization,  the  effective  electro-optic  coefficient  reff  in  our  measurement  is 
equal  to  r,3.  For  the  parameters  of  A.  =  633  nm,  0  =  1 1°,  n0  =  2.3,  m  =  1  and  d  =  1.5 
mm,  we  got  r13  =  48.5  pm/v  with  no  applied  field.  Then  with  an  applied  field,  the 
diffraction  efficiency  and  related  calculation  results  are  listed  in  Table  6-4.  We 
plotted  electro-optic  coefficient  versus  electric  field  and  fitted  the  data,  which  is 
shown  in  Fig.  6-6. 

The  fitting  results  gave  us 
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Table  6-4  Electric  Field  Dependence  of  the  Electro-Optic  Coefficient 
Eq  (  v/cm )  r|  ( xlO'3 )  Sin'1  ( rj 11 2 )  (  xl  0"2 )  Esc(  v/cm  )  reff 


r(  Eq)  /  r13  =  0.99308+  (  5.282  •  10  ~4  )  Eq  -  (  6.935  •  10  "7 )  Eq2  + 

(2.682  •  10  -I° )  E03  —  ( 3.289  •  10"14)E^4 
For  our  SBN:60  crystal,  the  effective  electro-optic  coefficient  Reff  is  given  by 
Reff  =  ne4  r33  cos20  -  n04  r13  sin20  for  extraordinary  rays  and  Reff  =  n04  r13  for  ordinary 
rays.  Experimentally  we  measured  the  gain  coefficient  of  the  weak  beam  with  an 
ordinary  polarization  and  an  extraordinary  polarization,  respectively.  Then  by 
calculating  the  ratio  of  two  gains,  we  can  get  the  ratio  of  r33  to  r13.  From  our 
measurement  the  ratio  of  r33  to  r13  is  3.21  for  our  crystal. 

Experimentally,  we  also  measured  the  electro-optic  coefficient  by  an  optical 
interference  experiment.  The  experimental  set  up  is  shown  in  Fig.  6-7. 


Fig.  6-7  Experimental  Setup  For  Measuring  The  Interference  Pattern 
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From  this  experiment,  the  electro-optic  coefficient  could  be  determined  by  measuring 
the  resultant  shift  of  the  interference  fringe  pattern  due  to  applying  an  electric  field 
across  the  crystal.  As  a  voltage  V  was  applied  to  a  crystal,  the  resultant  refractive 
index  change  in  the  crystal  would  be 

A11  =  (  r33ne3  '  ri3no3  )V/2d 

where  d  was  the  width  of  the  crystal.  This  refractive  index  change  caused  a 
difference  of  the  optical  path  between  two  beams,  which  was  equal  to  2*An*l  (1:  the 
length  of  the  crystal ).  Experimentally,  we  founded  that  the  fringe  pattern  moved  one 
space,  which  indicated  that  the  optical  path  difference  between  two  beams  was  equal 
to  A,  ,  when  the  voltage  across  the  crystal  was  625  v.  Based  on  this  data  we  can 
calculate  rI3  which  is  equal  to  38.1  pm/v  for  an  applied  voltage  of  625  v.  For  the 
same  voltage  of  625  v ,  from  the-data  in  table  6-4  we  got  r13(V0=625  v)  =  0.95*48.5  = 
46.1  pm/v. 

Comparing  these  two  values  of  r13  obtained  by  different  experiment  proved  that  our 
results  for  the  electric  field  dependence  of  the  electro-optic  coefficient  are 
dependable. 

Finally,  in  order  to  exam  the  accuracy  of  NA  value  which  we  obtained  by  using  the 
technique  of  determination  of  phase  shift.  We  also  did  an  experiment  to  measure  NA 
value  by  using  a  traditional  method  of  changing  the  angle  between  two  beams.  The 
basic  theory  and  measurement  steps  are  as  follows  (  Ref  [4] ). 

Due  to  two-beam  coupling,  the  weak  beam  Ij  experiences  an  increase  in  intensity 
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along  its  propagation  direction  which  is  (  see  Appendix  I ) 


Ilc  (l+m)eYLeff 
i!  l  +  meYLefr 


(6-6) 


where  m=  Iq,  /  I02,  a  ratio  of  the  weak  beam  intensity  to  the  strong  beam  intensity. 

Experimentally,  we  measured  Ilc,  I,  and  m,  and  then  calculated  y  Lcff. 

In  other  hand. 


Y  =R‘ 


27iRefrKBT  Kg 


XncosGe  l+(Kg/K0) 


(6-7) 


where  Reff  =  ne4  r33  cos20  -  n04  r13  sin20  for  an  extraordinary  ray,  Kg  is  grating 


2  4  7t  C  ^  N  pi-  ^ 

constant  and  Ko=  Kb”(“T")  =c(*fir)  provided  C  =  4n  eV  KBT,and  R  is  the 
measure  of  the  electro-hole  competition  in  the  formation  of  the  space-charge  field. 


n  -r 

By  defining  f(0)=cos0  /(cos20 - £*r-LJ-sin20)  we  will  have 

ne  f33 


K„ 


YLefrf(0)a 


"elT  - -  / 

C  /  £ 


By  plotting  the  data  of  y  Leff  f  ( 0 )  vs  Kg  and  fitting  the  data,  we  then  obtained  a  value 
of  Neff/e.  The  data  we  measured  and  calculated  are  listed  in  Table  6-5,  and  the 
fitting  curve  is  shown  in  Fig.6-8.  The  fitting  results  gave  us  the  value  of  Ncfl/e 


which  was 
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Fig.  6-8 

Neff/e  =  ( 4.59  ±  0.31 )  x  10 13  ( cm'3 ) 

Therefore,  from  the  data  we  have  obtained,  we  can  see  that 

1.  We  got  a  very  good  fitting  for  the  data  of  sin  tp  vs  E0  and  tan  tp  vs  E0  from  our 

phase  shift  measurement.  Also,  two  values  of  Nef(/€  obtained  in  two  different 
experiments  are  so  close.  All  these  provide  an  evidence  that  the  new  technique  of 
measuring  NAdeveloped  by  us  is  accurate  and  dependable. 

2.  Based  on  the  experimental  data,  we  got  fitting  curves  for  the  electric  field 
dependance  of  both  the  dielectric  constant  and  the  electro-optic  coefficient  which 
were  shown  in  Fig  6-5  and  6-6.  In  Ref.  [5]  the  author  also  plotted  these  two  curves 
of  €  vs  E0  and  r  vs  Eq.  Ours  have  the  same  shape  with  theirs.  This  varifies 
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indirectly  the  validity  of  our  new  developed  technique. 


IV.  Calculations  Of  Electron-Hole  Competition 


If  electrons  and  holes  occur  simultaneously  in  the  charge  transport  process, 
there  will  be  a  electron-hole  competition  in  the  formation  of  the  space-charge  field 
which  will  lead  to  a  reduction  in  strength  of  the  photorefractive  effect.  Suppose  both 
electrons  and  holes  can  be  excited  from  (  and  to  )  a  single  impurity  band,  then  by 
employing  band-transport  model  in  Ch.  3,  a  first-order  differential  equation  for  the 

time  evolution  of  the  space-charge  field  can  be  derived  which  is  (  see  Ref.  [  6  ] ) 

dE,(t)  D  , ,  ,  nmepR[l+k(k-iV)/kn2]  nhneph [  1  +k(k+i  V)/kn2] 

'  dt_’  =  E'(t)  e  [l+k(k-iV)/Ke2]  +  e  [  1 +k(k+i  VVK,,2]  J 

,  imknT  f  n^ep.fk-iV)  nhneph(k+iV)  , 

+  e  1  6  [l+k(k-iV)/K  2]  '  e  [l+k(k+iV)/K  2] ;  1  0-6 ' 


where  ne0  =  scI0(ND-NA)/yeNA  and  nh0  =  shI0NA/yh(ND-NA)  are  zeroth-order  free-carrier 
density. 

Ke-2  =  kBTpe/eyeNA,  Kh-2  =  k^/ey^-NJ,  kQ2  =  e2NA(ND-NA)/e  kBTND  and 
V  =  eE(/kBT.  Ke_1  and  Kh'‘  are  average  distance  travelled  by  electrons  and  holes 
between  excitation  and  recombination,  Eq  is  the  externally  applied  field. 

If  the  charge  carrier  is  only  electron,  then  p,=0,  Kh‘‘=0,  and  K,,— ><».  If  the  charge 
carrier  is  only  hole,  then  pe=0,  Ke-’=0,  and  Ke— >°o. 

If  we  write  eq.  (6-8)  as 
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-dEj  (t)/dt  =  AE ,  (t)  +  B  (  6-9  ) 

The  solution  for  eq.  (6-16)  is 
E1(t)  =  (Be-Al-B)/A 

At  equilibrium,  the  steady-state  space-charge  field  is 

esc=-b/a 


As  no  voltage  is  applied  to  a  crystal,  v  is 

A=(efe)(l+lttM(.  i;g/^- 

B  -(imknT  ,  neope 

l“^^)(e/fe)k(  1+kVK* 


equal  to  zero  in  eq.  (6-8). 
,  ^hol^h _ \ 

l+k2/Kh2  J 

-  ^hol^h  s 

l+k2/Kh2  ' 


Then  we  have 


Thus,  the  steady-state  space-charge  field  is 


^sc  E  V 


-K- 


^co.^c 


^hoM-h 


/(- 


^eol^e 


l+k2/Ke2 


l+k2/K02  v  l+k2/Ke2 
_ ,  ^hoM-h  \ 


l+k2/Kh2 


l+k2/Kh2 


or 


p  _  imknT  k 

-L'co - -  A11  ’  — . . 


*  l+k2/k02  R 


where  R  =  ( 1-C  )/( 1+C ) 

with  C  =  shNA(kg2+Ke2)/[sc(ND-NA)(kg2+Kh2)] 

If  the  charge  carrier  is  only  electron  C=0,  R=l,  and  if  the  charge  carrier  is  only  hole 
C— R=-l .  Therefore,  R=1  or  -1  indicates  no  electron-hole  competition.  From 
our  data  in  table  6-5  we  plotted  Kg/y  vs  Kg2and  fitted  the  points.  The  results  gave  us 
a  straight  line  which  is  shown  in  Fig.  6-9.  The  slope  and  the  y-intercept  of  the 
line  enable  us  to  determine  the  electro-hole  competition  R,  which  is  found  to  be 
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1.4 


kg2  (xIO4  1/cm2) 

Fig.  6-9  Determination  of  Electron-Hole  Competition  Factor 


0.98.  This  indicates  that  there  is  no  electron-hole  competition  in  the  formation  of  the 
space-charge  field. 


V.  Conclusion 

We  have  demonstrated  a  new  technique  for  the  measurement  of  the  trap  density 
in  photorefractive  materials.  This  technique  has  been  shown  to  give  the  same  result 
for  the  trap  density  as  the  conventional  two-beam  coupling  technique  but  is  more 
convenient  to  carry  out  in  practice.  The  technique  presented  here  also  yields 
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the  electric  field  dependence  of  both  the  dielectric  constant  and  the  electro-optic 
coefficient.  In  addition,  an  electro-hole  competition  was  determined  based  on  our 
data  of  gain  versus  grating  constant.  The  result  indicates  that  there  is  no  electron-hole 
competition  for  our  crystal. 


VII.  GRATING  WRITING  AND  ERASURE  TIMES 


In  developing  an  effective  photorefractive  optical  limiter  it  is  important  to  have 
in  hand  a  good  understanding  of  the  physics  of  the  response  time  of  the  material. 
Therefore,  part  of  our  research  effort  has  been  focused  on  developing  a  clear  picture  of 
the  writing  and  erasing  photorefractive  response  times.  Our  results  are  in  strictly 
excellent  agreement  with  our  photorefractive  model  although  some  are  intuitively 
surprising.  For  example,  we  have  found  the  I'1  dependence  to  be  vigorously  followed 
for  our  SBN  samples  while  we  have  found  that,  contrary  to  suggestions  in  the 
literature,  the  writing  and  erasing  times  are  not  equal  for  equal  incident  intensities. 


Investigation  On  Grating  Formation  Time  And  Erasure  Time 
Analytical  Solution  For  Grating  Erasure  Time 

When  a  grating  is  being  written  the  total  current  due  to  the  electrons  in  the 
conduction  band  can  be  expressed  as  j  =  j0+j1eikx,  where  j0=e(in0E0  depends  on  the 
applied  field  and  j,  has  two  parts  jD  and  jE.  The  jD  part  is  a  current  due  to  the  diffusion 
and  the  jE  is  a  current  due  to  the  formation  of  the  space-charge  field.  In  this  physical 
picture,  the  direction  of  net  electron  flow  is  the  same  as  the  electric  field  direction.  But 


Electron  Distribution 
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when  a  grating  is  erasing  the  situation  is  different  from  the  writing  case. 

1 .  The  current  still  has  two  parts  j0  and  j,,  however,  the  direction  of  net  electron 
flow  is  now  opposite  to  the  direction  of  the  space-charge  field.  Therefore,  we 
have  to  write  the  current  as  j  =  j0  -  j  je,kx  where  the  minus  sign  can  be  explained  as 
a  1 80°  phase  shift. 

2.  In  the  grating  erasing  case  the  electron  are  driven  by  the  space-charge  field.  The 
part  of  j,  due  to  diffusion  can  often  be  neglected. 

3.  During  the  period  of  erasing  a  grating  the  optical  intensity  is  uniform,  and  can  be 
written  as  I  =  I0. 

Based  on  the  above  consideration  the  four  basic  equations  describing  the 

photorefractive  effect  in  a  grating  erasing  case  should  be  written  as 


iin  _  dNp  __  \ 

3t  9t  -  e  5x 

3N* 

— D=SI(N-ND)-YRnND 
j  =  ep.nE 

M=4ffi(n  +  NA_N5 

Then  the  assumptions  are 


(7-1) 

(7-2) 

(7-3) 

(7-4) 


I  =  Io 


(7-5) 


n  =  n0+  nj(t)  eikx 
ND=N>N^(t)eikx 

E  =  E0+  E,(  t )  eikx 

j=Jo-ji(t)eikx 


(7-6) 

(7-7) 

(7-8) 

(7-9) 


By  substituting  ( 7-5  )  -  ( 7-9  )  into  ( 7-1 )  -  ( 7-4 )  we  can  obtained  a  second-order 
differential  equation  which  is 


a2Ej  f  _1_  _  J_  _  _L )  .  / _ L_ 

9t2  ^di  Xg  0t  ^di^R 


i 

X,XE 


•)  E  =  0 


(7-10) 


The  solution  for  ( 7-10 )  is 


E,(  t )  =  E^e  '^e.  e'“‘  where 


is  the  grating  erasure  time. 


(7-11) 
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In  the  writing  case,  the  grating  writing  time  (  or  say  "  grating  growing  time"  )  is 
given  by  eq.  (4-5),  which  is 


In  the  case  that  ^<<5^  eq.  (4-5)  can  be  written  as 


(7-12) 


If  we  now  compare  (  7-11  )  with  (  7-12  ),  it  is  clear  that  the  grating  erasure  time  is 
longer  than  the  writing  time.  This  is  true  even  in  the  case  of  zero  field  but  is 

particularly  important  for  the  applied  field  case  since  xE  can  be  on  the  order  of  xR 
and  xdi  can  be  10  times  longer  than  x,. 


Experimental  Measurement  For  Grating  Writing  Time  And  Erasing  Time 

Experimentally,  we  measured  the  writing  time  and  the  erasing  time  of 
photorefractve  index  gratings  for  SBN:60  with  Cerium  as  a  dopant.  The 
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experimental  set  up  is  shown  in  Fig.  7-1. 


Determination  for  the  varying  function  of  the  grating  erase  with  time 
asingle  exponential  decay  function 
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In  order  to  measure  the  grating  erasing  time  we  must  first  know  the  function 
of  how  the  grating  erasure  varying  with  time.  We  measured  a  diffraction  signal  for  a 
SBN  sample  as  a  function  of  time.  The  results  are  shown  in  Fig.  7-2,  which 
indicated  that  the  grating  erase  is  a  single  exponential  time  decay  function. 
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Fig.  7-2  Single  Exponential  Decay  Function 


Measurement  for  grating  growing  time  and  erasing  time 
Our  experimental  parameters  were:  wavelength  =  458  nm,  total  laser  power  =  25 
mw,  the  full  crossing  angle  between  two  beams  =  10°,  and  the  intensity  ratio  of  two 
beams  =  1:100.  Experimentally,  the  weak  beam  was  detected.  At  first  the  energy 
exchange  signal  was  traced  by  DSA  oscilloscope,  from  which  a  grating  growing  time 
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xg  can  be  measured.  After  the  refractive  index  grating  being  formed  inside  the 
crystal  we  blocked  the  detected  beam  and  then  the  diffraction  signal  due  to  the  strong 
beam  was  traced  by  DSA  oscilloscope,  from  which  a  grating  erasure  time  xe  can  be 
measured. 

The  curve  of  grating  formation  is  shown  in  Fig.  7-3,  and  that  of  grating  erase  is  in 
Fig.  7-4.  The  measurement  results  for  the  grating  growing  time  and  the  erasing  time 
are  xg=2.66  s  and  xe=4.22  s,  respectively,  which  showed  that  the  grating  erasure  time 
is  longer  than  the  grating  formation  time.  Quantitatively,  in  this  measurement  the 


Xg=  2.66  s 
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Fig.  7-4  Grating  Erasure  Time 

grating  erasure  time  is  about  1.6  times  of  the  grating  formation  time. 

Another  measurement  of  the  grating  writing  time  x  and  the  erasing  time  xe  was 
performed  for  BaTi03  crystal  with  different  kg  values.  The  measurement  results  are 
(1)  \!\  =  1.5  for  crossing  angle  of  14°;  (2)  xe/xg  =  1.9  for  crossing  angle  of  40°. 
These  results  agree  with  the  theoretical  predictions  that  (1)  the  grating  erasure  time  is 
longer  than  the  grating  growth  time  and  (2)  the  ratio  of  the  grating  erasure  time  to  the 
grating  growth  time  is  increased  with  kg  increased. 

Comparison  Of  Experimental  Data  With  Theoretical  Calculation 

In  the  case  of  no  applied  field,  eq.  (7-1 1)  gave  an  erasing  time  of  xc  =  xdjand  eq. 
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(7-12)  gave  a  growth  time  of 


1  +(-*)(-f) 


For  the  parameters  of  X=458  nm,  29=10°,  NA=D1016  cm'3  and  ND=10NA,  we  can 
calculate  the  ratio  of  grating  erasure  time  to  formation  time  which  is  1.8.  This 
theoretical  value  is  close  to  the  experimental  result  for  the  ratio  of  1.6  which  we 
measured.  Therefore,  the  experimental  results  conform  the  theoretical  analysis  for 
the  photorefractive  grating  erasure  time  which  we  derived  in  the  beginning  of  this 
chapter.  Both  theoretical  and  experimental  results  do  indicate  that  the  grating 
erasure  time  is  longer  than  the  growth  time. 


From  eq.  (7-11)  and  (7-12)  we  can  obtain 

_L  = - Ln  -  xi  lF  ~i 

T«  1-  (^X-^-) 

te  ti 


=  (1 +Akg  +  Bkg2)/(1 -Bkg2)  (7-13) 

where  A  =  pE^.  /  x,  and  B  =  \i%\\  /  x,. 

From  eq.  (7-13),  we  can  see  that  the  ratio  of  the  grating  erasure  time  to  the 
growth  time  is  increased  with  kg  increased. 


Dependence  Of  Grating  Erasure  Time  On  Applied  DC  Filed 


In  the  last  section  we  have  seen  that  in  the  case  of  no  applied  field  the  grating 


erasure  time  is  longer  than  the  formation  time,  both  theoretically  and  experimentally. 
We  have  measured  the  grating  erasure  time  without  field,  which  was  4.22  s  (  see 
Fig.4-7  ).  But  what  happens  when  a  DC  field  is  applied?  According  to  our 
theoretical  analysis  given  by  eq.(7-ll),  the  grating  erasing  time  should  get  longer 
with  a  DC  field  than  that  with  no  field.  In  order  to  verify  the  theoretical  prediction, 
experimentally,  we  measured  a  grating  erasure  time  with  a  DC  field  applied  to  the 
crystal.  The  experimental  parameters  were  same  as  those  described  in  last  section. 
When  a  DC  field  of  6  kv/cm  was  applied  to  the  crystal,  a  grating  erasure  time  of  6.4  s 
was  measured.  This  erasing  time  of  6.4  s  is  obviously  longer  than  the  erasing  time  of 
4.22  s  with  no  field  obtained  in  last  section  (  Fig.4-7  ). 

Then  with  eq.(  7-11 ),  we  calculated  the  theoretical  value  of  grating  erasing  time 
with  a  DC  field  of  6  kv/cm.  The  calculated  result  for  grating  erasure  time  with  a  DC 
field  of  6  kv/cm  was  7.0  s.  This  is  close  to  the  experimental  result  of  6.4  s. 
Therefore,  the  theoretical  prediction  for  the  grating  erasure  time  given  by  eq.(  7-11) 
is  dependable.  Both  theoretical  and  experimental  results  indicates  that  the  grating 
erasure  time  depends  on  the  applied  field.  In  the  case  of  a  DC  field  applied  to  a 
crystal  the  grating  erasure  time  is  getting  longer  than  that  in  the  case  of  no  field. 

Dependence  Of  Grating  Erasure  Rate  On  Optical  Intensity 

Photorefractive  crystals  have  been  used  to  demonstrate  several  novel  operations 


118 


in  real-time  optical  data  processing.  Understanding  the  time  dependent  behavior  in 
these  materials  is  important  both  for  predicting  the  response  of  the  photorefractive 
crystals  and  for  determining  how  to  modify  them  to  improve  response  (  Ref.  [2]-[5] 
)• 

In  this  chapter  we  will  present  the  data  which  show  the  dependance  of  light- 
induced  erasure  rate  of  a  refractive  index  grating  on  optical  intensity.  Our 
measurement  was  performed  for  SBN  crystals  and  BaTiOs  crystals,  respectively. 
Different  behavior  of  the  dependance  were  obtained  for  these  two  types  of  crystals. 
The  curve  fitting  results  showed  that  in  SBN  the  speed  of  erasure  of  a  photorefractive 
index  change  is  linearly  dependent  on  the  intensity  of  the  erasing  beam  with  and 
without  an  applied  field.  However,  in  BaTi03  this  erasure  rate  of  a  photorefractive 
grating  scaled  sublinearly  with  optical  intensity,  with  and  without  a  field. 

Experimental  Setup  And  Data  Collection 

The  experimental  setup  is  shown  in  Fig.  7-1.  Two  writing  beams  wrote  a 
photorefractive  grating  inside  a  crystal.  The  third  beam  which  bisected  two  writing 
beams  was  used  to  erase  the  grating.  It  is  called  erasing  beam.  Experimentally,  we 
detected  the  weak  beam.  After  the  grating  being  formed  inside  the  crystal  we 
blocked  the  detected  beam,  and  at  the  same  time  turned  on  the  erasing  beam.  Then  a 
diffraction  signal  from  another  beam  was  sent  to  PMT  and  measured  by  DSA.  By 
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changing  the  intensity  of  erasing  beam  we  are  able  to  measure  the  dependance  of  the 
grating  erasure  rate  on  the  optical  intensity.  Our  measurement  were  performed  for 
two  SBN  thin  samples,  two  BaTi03  thin  samples  and  one  BaTi03  thick  sample.  The 
experimental  data  of  grating  erasure  time  versus  erasing  beam  intensity  are  listed  in 
Table  7-1  to  Table  7-3.  The  fitting  results  are  shown  in  Fig.  7-5  to  Fig.  7-11.  From 
the  fitting  curves  we  found  that  there  is  a  very  good  linear  dependance  of  the  grating 
erasure  time  on  optical  intensity  for  SBN  samples.  But  for  BaTi03  samples  the 
behavior  shown  are  apparently  sublinear.  The  fitting  results  gave  an  r0'67 
dependence  for  one  BaTi03  thin  sample,  f0  72  dependence  for  another  BaTi03  thin 
sample  and  r(0  8  0,84)  dependence  for  the  BaTi03  thick  sample. 

Table  7-1  Data  of  Measuring  the  Linear  Behavior  for  SBN  Crystals 
Thin  Sample  #1 

Neutral  density  filter  Ierase(mw)  x’(s)  Tdark(s)  l/x=l/x’- l/xdark  ( i/s) 
in  erasing  beam 


0 

75 

0.052 

790 

19.23 

1 

5.5 

0.67 

790 

1.491 

2 

0.55 

6.7 

790 

0.148 

3 

0.055 

70 

790 

0.013 

4 

0.0055 

376 

790 

0.00139 
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Thin  Sample  #2 


Neutral  density  filter  Ierase(mw)  x’(s)  x^s)  l/x=l/x’-l/xdark  (l/s) 
in  erasing  beam 


0 

20 

0.1503 

75.633 

6.6419 

1 

2 

1.52 

75.633 

0.6447 

2 

0.2 

13.47 

75.633 

0.061 

3 

0.02 

50.706 

75.633 

0.0065 

4 

0.002 

71.369 

75.633 

0.00079 

0.001  0.01  o.l  l  10  100 

Optical  power  in  erasing  beam  (  mw  ) 


Fig.  7-5 


Fig.  7-6 

Table  7-2  Data  of  Measuring  the  Linear  Behavior  for  BaTi03  Crystals 
Thin  Sample  #1 

Wmw>  T’(s)  'W8)  tn/x=i/t’-i/tdark(i/S)] 


0.8 

19.0 

40 

36.2 

2 

13.4 

40 

20.2 

3.2 

10.8 

40 

14.8 

8 

6.50 

40 

7.76 

32 

2.92 

40 

3.15 

80 

1.62 

40 

1.69 
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Thin  Sample  #2 

Wt  mw  )  T’  (S)  Tdark  (  s  )  x[l/T=l/t’-l/xdark(l/s)] 


BaTiOo  Thin  Sample  2  ~4 - I  °'724 

J  erase 


Optical  power  in  erasing  beam  (  mw  ) 
Fig.  7-8 


Table  7-3  Data  of  Measuring  the  Linear  Behavior  for  BaTi03  Thick  Sample 
BaTi03  Thick  Sample  ( #3  ) 

Table  A:  Full  crossing  angle  20=14°  Table  B:  Full  crossing  angle  20=40° 
Intensity  ( mw )  l/xerase  (  1/s  )  Intensity  ( mw  )  l/xerase  ( 1/s ) 


217 

0.6235 

200 

0.5665 

138.23 

0.3715 

100 

0.254 

84.41 

0.221 

78 

0.194 

42.75 

0.113 

39.4 

0.1125 

33.2 

0.092 

30.6 

0.073 

13.67 

0.0615 

19.8 

0.048 

8.46 

0.0355 

7.8 

0.021 

4.25 

0.0205 

3.92 

0.016 

3.28 

0.018 

3.02 

0.014 
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Optical  power  in  erasing  beam  ( mw 

Fig.  7-9 


BaTi03  sample  3 

Full  crossing  angle:  40  degree 


Optical  power  in  erasing  beam  (  mw  ) 


Table  7-3  Data  of  Measuring  the  Linear  Behavior  for  BaTi03  Thick  Sample 
BaTi03  Thick  Sample  ( #3  )  Table  C:  Full  crossing  angle  20=80° 
Optical  Intensity  (  mw  )  1/Terase  (  1/s  ) 


Explanation  For  Sublihear  Photorefractive  Response  Time  in  BaTi03 

From  the  band  transport  model  the  photorefractive  response  time  x  is  governed 
by  eq.  ( 4-5  ),  which  gives  x  is  proportional  to  xdi,  the  dielectric  relaxation  time.  This 
xdi  is  inversely  proportional  to  n0,  the  zeroth-order  electron  number  density.  In  the 
single-donor-trap  model  n0  is  the  only  intensity-dependent  parameter  that  affects  eq. 
(  4-5  ).  It  is  simply  proportional  to  the  laser  intensity,  and  consequently  the 
response  time  x  is  expected  to  vary  inversely  with  intensity.  The  results  of  our 
measurement  for  the  dependance  of  the  grating  erasure  rate  on  optical  intensity 
showed  a  good  agreement  with  the  theoretical  prediction.  But  this  linear  behavior 
only  happened  to  crystals  the  donor  number  density  of  which  greatly  exceeds  that  of 
accepters  (  ND  »  NA  ).  For  .this  type  of  crystals  which  were  defined  as  type- A 
crystals  (  Ref.  [2]-[5]  ),  the  previous  models  of  the  photorefractive  effect  assumed 
that  either  the  donors  or  the  accepters  (  but  not  both  )  took  part  in  charge  transport, 
and  neglected  thermal  excitations. 

However,  for  another  type  of  crystals  which  have  ND  <  NA  and  are  defined  as 
type-B  crystals  (  Ref.  [2]- [5]  ),  the  band  transport  model  should  be  extended  to 
include  secondary  photorefractive  centers.  These  are  trapping  centers  that  are  highly 
ionized  at  room  temperature  but  become  populated  by  acquiring  free  charge  carriers 
generated  by  the  interaction  of  the  laser  light  with  the  deep-level  impurities.  The 
important  feature  of  this  model  is  the  finite  thermal  ionization  rate  of  the  secondary 
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centers,  which  influences  the  free-carrier  and  occupied  trap  concentrations. 
Theoretically,  one  more  rate  equation  for  accepters  should  be  added  to  the  original  four 
basic  equations  (3-1)  to  (3-4),  then  the  nonlinear  behavior  between  the  grating  erasure 
rate  tne  the  optical  intensity  for  BaTi03  could  be  explained  (Ref.  [6]). 

Conclusionsn 

Our  results,  both  theory  and  experiment,  show  that  the  writing  and  erasure  time 
for  photorefractive  materials  are  not  the  same.  However,  they  are  in  practice 
sufficiently  close  considered  to  be  on  the  same  order.  At  the  same  time  an  examination 
of  response  time  versus  intensity  dramatically  supports  a  I1  dependence  for  SBN.  As  a 
result  our  assumption  for  a  photorefractive  limiter  based  on  SBN  is  validated. 
However,  this  assumption  is  not  true  for  all  materials. 
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VIII.  SELF-FOCUSING  AND  DEFOCUSING 

One  of  the  most  dramatic  photorefractive  phenomena  observed  in  the  laboratory 
has  been  self-focusing  and  self-defocusing.  Either  phenomena  can  be  used  together 
with  an  imaging  lens  and  aperture  as  a  limiter. 

Fig.  8-1  shows  a  top  view  of  a  crystal  with  a  10  micron  diameter  beam  passing 
through  it.  Without  an  applied  voltage  the  beam  is  clearly  seen  to  diverge  due  to  normal 
diffraction.  With  a  small  applied  voltage,  however,  the  beam  is  seen  to  defy 
diffraction.  This  effect  is  due  to  the  self-focusing  photorefractive  phenomena. 

As  ,in  the  case  of  Kerr  focusing  and  defocusing,  photorefractive  focusing  and 
defocusing  can  also  be  physically  understood.  To  illustrate  the  photorefractive  effect, 
we  can  first  consider  two  plane  waves  of  light  overlapping  in  a  crystal  producing  an 
optical  interference  pattern  as  shown  in  Figure  8-2.  In  the  bright  regions  of  the 
interference  pattern  carriers  are  excited  into  the  conduction  band.  The  excited  carriers 
then  diffuse  or  drift  and  are  finally  trapped  in  the  dark  regions  of  t  interference  pattern. 
The  resulting  charge  separation  in  turn  generates  a  space  charge  electric  field.  The  drift 
and  diffusion  process  continues  until  equilibrium  is  reached  where  the  diffusion  or  drift 
current  is  exactly  balanced  by  the  current  generated  by  the  induced  space  charge  field. 
In  this  way  the  magnitude  of  the  field  is  simply  determined  by  the  value  necessary  to 
balance  the  diffusion  or  drift  current.  The  space  charge  field  can  then  distort  the  lattice 
and  produce  via  the  electro-optic  an  index  change  given  by 

n3 

An=~  rc(rEsc  (8-1) 

ere  refr  is  the  effective  electro-optic  coefficient  of  the  material,  Esc  is  the  induced  space 
charge  field,  and  n  is  the  unperturbed  index  of  refraction  of  the  material. 

The  resulting  induced  index  change  can  then  be  used  in  Maxwell’s  equation  to 
predict  the  propagation  behavior  of  the  two  overlapping  laser  beams.  The  result  is  that 
the  induced  index  causes  a  coupling  between  the  two  beams  which  can  be  written  as 
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(8-2) 


form  a  physical  picture  of  diffraction  that  is  based  on  “watching”  each  Fourier 
component  propagating  through  the  material.  Since  each  component  has  a  different  k- 
vector  projection  direction,  the  relative  phase  between  Fourier  components  changes  as  a 
function  of  propagation  distance  z.  Consequently,  the  sum  of  the  Fourier  components 
produces  a  different  wave  form  at  each  z  or  propagation  position.  In  the 
photorefractive  picture  each  Fourier  component  produces  an  interference  pattern  with 
each  and  every  other  Fourier  component.  The  result  is  that  each  component,  therefore, 
“sees”  a  modified  index  of  refraction  which  is  determined  by  summing  the  index 
modification  produced  between  a  given  Fourier  component  and  every  other  component. 
When  the  low  frequency  Fourier  components  “see”  a  lower  index  than  the  higher 
frequency  components  due  to  the  coupling,  focusing  is  induced.  Figure  8-3(b) 
Likewise,  when  the  low  frequency  components  “see””  a  higher  index  than  the  higher 
frequency  components,  defocusing  is  induced  (Figure  8-3(c).  In  practice,  the  sign  of 
the  index  change  depends  on  the  sign  of  an  applied  external  field  so  that  focusing  is 
induced  with  an  applied  field  along  the  c-axis  direction  while  defocused  field  opposite 
the  c-axis  direction.  On  interesting  possibility  not  discussed  here  occurs  when 
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diffraction  is  exactly  compensated  by  photorefractive  focusing.  In  this  case,  shape 
pursuing  propagation  or  soliton  formation  is  observed.  This  picture  neglects  energy 
coupling  which  would  result  in  an  amplitude  change  of  the  Fourier  components  and 
would  be  important  for  a  more  complete  analysis.  In  our  work  the  effect  of  energy- 
coupling  was  minimized  by  the  choice  of  input  beam  diameter,  the  value  of  applied 
field,  and  the  time  at  which  the  measurement  was  made. 

The  basic  apparatus  consisted  of  a  cw  argon-ion  laser  and  a  5  mm  x  5  mm  x  6 
mm  strontium  barium  niobate  (SBN)  crystal  with  0.01%  by  weight  rhodium  dopant. 
The  cw  argon-ion  laser  wavelength  was  457  nm  and  its  output  beam  diameter  was  1.5 
mm.  A  schematic  diagram  of  the  apparatus  is  shown  in  Figure  8-4.  The  output  beam 
was  directed  onto  a  10  cm  focal  length  lens  and  the  SBN  crystal  was  placed  2.6  mm 
beyond  the  beam  waist  with  2  w0  =  33  pm.  The  beam  diameter  at  the  SBN  crystal 
entrance  face  was  75  pm.  The  crystal  was  oriented  with  its  c-axis  in  the  horizontal 
plane  and  perpendicular  to  the  propagation  direction  of  the  incoming  laser  light.  The 
polarization  of  the  incoming  light  could  be  varied  using  a  polarization  rotator  but  was 
initially  chosen  to  be  along  the  c-axis  (extraordinary  polarization).  The  beam  diameter 
throughout  the  crystal  was  measured  using  an  imaging  system  consisting  of  an  imaging 
lens  and  a  two-dimensional  detector  array.  The  input  face  of  the  6  mm  long  SBN 
crystal  was  well  beyond  the  Rayleigh  range  of  1.5  mm  from  the  beam  waist  formed  by 
the  focusing  lens.  The  imaging  system,  therefore,  imaged  the  beam  spot  at  the  SBN 
entrance  face  with  some  magnification  onto  the  detector  array.  As  the  imaging  lens  and 
the  detector  array  are  moved  away  from  the  SBN  crystal,  different  cross  sections  of  the 
Gaussian  beam  are  then  imaged  onto  the  array.  In  this  manner,  the  beam  diameter  at 
different  locations  throughout  the  SBN  crystal  was  monitored.  The  magnification  of 
the  imaging  system  was  determined  by  placing  a  thin  aperture  on  the  crystal  exit  (and 
entrance)  face  and  imaging  the  aperture  onto  the  detector  array.  Using  the  known  value 
of  the  reference  aperture,  the  magnification  was  determined  to  be  about  15.6  and  the 
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positions  of  the  exit  and  entrance  faces  of  the  SBN  crystal  were  located.  Using  this 
information  the  horizontal  cross  section  of  the  incident  beam  on  the  entrance  and  exit 
faces  of  the  crystal  was  determined.  Beam  diameter  and  viergence  data  were  taken  at  a 
time  when  these  quantities  reached  a  maximum  value. 

Figure  8-5  shows  the  effect  of  applied  voltage  on  the  exiting  beam  diameter 
while  Figure  8-6  shows  its  effect  on  the  divergence  of  the  incident  beam.  As  shown  in 
these  two  figures  the  photorefractive  defocusing  is  a  dramatic  effect.  Similar  results 
have  been  observed  for  focusing.  Only  the  sign  of  the  applied  d.c.  electric  field  was 
reversed  when  producing  focusing  or  defocusing  effects.  Defocusing  effects  are  seen 
when  the  field  is  opposite  the  c-axis  direction. 

Conclusions 

Our  results  clearly  demonstrate  the  potential  of  photorefractive  self-focusing  and 
defocusing  as  an  optical  limiter.  Both  are  seen  to  be  large  and  dramatic  effects  to  occur 
with  little  or  low  laser  power.  Although  we  have  not  had  the  opportunity  to  study  this 
effect  e  have  determined  the  magnitude  and  speed  of  the  effect  and  its  independence  on 
intensity.  That  is,  with  the  same  applied  voltage  the  degree  of  focusing  and  defocusing 
is  observed  to  be  independent  of  intensity  over  four  orders  of  magnitude,  as  expected. 
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CONTROL  OF  NA  AND  ND 

As  has  been  pointed  on  numerous  occasions  in  this  report,  control  of  NA  and  ND  is 
important  to  the  design  of  a  photoreffactive  optical  limiter.  Variations  in  the  response  time 
and  the  O.D.  and  NA  and  ND  are  large.  One  potential  technique  for  control  of  both  of 
these  variables  is  to  use  separate  dopants  in  an  attempt  to  control  them  independently.  For 
this  reason  we  have  investigated  the  behavior  of  double-doped  crystals.  In  particular,  we 
have  focused  on  the  Cr-Mo  combination.  Four  crystals  were  studied  and  the  results  are 
shown  below. 

Cr,  Mo  DATA 

(a)  .01  Cr  and  .005  Mo 

a  L -.66 

r  ~  1.68  sec 
OD  -  1.6 

all  data  at  1W/CM2  at4880A 

(b)  .005  Cr  and  .01  Mo 

a  L  -.33 

r  —  1.68  sec 
OD  -  1.3 

all  data  at  1 W/CM2  at  4880A 

(c)  .005  Cr  and  .005  Mo 

a  L  -.33 

r  -  .8  sec. 

OD  -  1.2 

all  data  at  1  W/CM2  at  4880A 
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(d)  .01  Cr  and  .01  Mo 

a  L  -.55  (.65  for  second  cube) 

t  ~  1.7  sec. 

OD  ~  .3 

all  data  at  1 W/CM2  at  4880A 

The  data  from  crystals  (a),  (b),  and  (c)  supports  the  idea  that  only  Cr  controls  the 
crystal  absorption.  That  is  Cr  plays  the  role  of  ND.  For  crystals  (a)  and  (b)  the  Cr 
concentration  changed  by  a  factor  of  2.  Meanwhile  the  total  dopant  concentration  remained 
constant.  This  indicates  that  Cr  alone  plays  the  role  of  the  optical  absorber.  Likewise 

crystal  (c)  has  the  same  Cr  concentration  as  crystal  (b).  For  these  two  crystals  aL  or  the 

absorption  was  determined  to  be  the  same.  These  three  crystals  therefore  offer  strong 
support  that  Cr  alone  plays  the  role  of  ND.  Crystal  (d)  behaves  more  like  a  Cr 
concentration  of  .0075  and  a  Mo  concentration  of  .0075.  Perhaps  this  was  tcut  from  a 
different  boule  location. 

The  same  crystals  also  indicate  that  both  Cr  and  Mo  are  playing  the  role  of  NA. 

That  is,  the  response  time  for  both  crystals  (a)  and  (b)  measure  to  be  the  same.  This  would 
be  the  case  if  the  trap  density  were  also  the  same  in  both  crystals.  Meanwhile  crystal  (c) 
has  a  clear  faster  response  time.  The  total  dopant  concentration  is  less  in  this  case  than  for 
crystal  (a)  or  (b).  Likewise  the  OD  for  crystal  (c)  is  less  than  that  for  crystals  (a)  and  (b) 
supporting  a  smaller  NA  for  crystal  (c).  The  fourth  crystal  looks  like  a  Cr,  Mo 
concentration  of  .0075  each  or  a  total  again  of  .015. 

The  data  from  crystal  (a),  (b),  (c),  and  (d)  therefore  support  the  idea  that 

*  Cr  plays  the  role  of  ND 

•  Cr  and  Mo  play  the  role  of  NA 


141 


In  addition  to  the  double-doped  samples  we  also  examined  the  two  types  of  dopants 
separately.  For  Mo  alone  we  found  for  four  samples: 

(a)  .005  Mo 

a  1  -  3.2 

x  -  1.2  sec 
OD  -  1.0 

All  data  at  lW/cm2  at  4880A 

(b)  .01  Mo 

al~  1.7 

t  ~  4  sec 
OD~  1.6 

All  data  at  1  W/cm2  at  4880A 

(c)  .015 

a  1  -2.9 

r  -  .228  sec 
OD  -  2.6 

All  data  at  lW/cm2 

(d)  .015 

a  1  -2.3 

t  -  1.26  sec 
OD  -  2.3 


All  data  at  lW/cm2  at  4880A 


(e)  .025 

a  1  -1.7 

t  -  1.98  sec. 

OD-3 

All  data  at  lW/cm2  at4880A 
(f)  .025 

al  -1.4 

t  -  2.35  sec. 

OD-3 

All  data  at  lW/cm2  at  4880A 
While  for  Cr  alone  we  found 

(a)  .005  Cr 

a  L  -  .35 

t  -  2.6  sec 
OD-  1 

all  data  at  1  W/cm2  at  4880A 

(b)  .01  Cr 

a  L  -  .69 

r  -  2.6  sec 
OD  -  1.3 

all  data  at  1  W/cm2  at  4880A 

(c)  .015  Cr 

a  L  -  .73 
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r  ~  .390  sec 


OD-1 

all  data  at  1  W/cm2  at  4880A 

(d)  .025  Cr 

a  L  ~1.2 

t  ~  .515  sec 
OD~  1 

all  data  at  1  W/cm2  at  4880A 

(e)  .015  Cr 

a  L -.8 

t  ~  .375  sec. 

OD-  1.3 

(f)  .025  Cr 

a  L~2 

r  ~  .455  sec. 

OD-1 

That  is,  for  Mo  we  have 

Data  from  crystal  (a),  (b),  (c),  and  (d)  suggest  that  Mo  plays  the  role  of  NA.  By 
increasing  the  concentration  of  Mo  we  increase  the  response  time  indicating  that  NA  is 

increasing.  Meanwhile  the  aL  is  decreasing!.  This  would  suggest  that  the  higher 

concentration  of  Mo  is  lowering  the  ND  allowed  in  the  crystal.  The  higher  NA  value  is  also 
supported  by  a  higher  OD  value  for  crystals  with  higher  Mo  concentration.  It  is  also 

interesting  that  the  ratio  of  aL  for  the  two  .015  crystals  is  2312.9  or  .80  while  the  ratio  of 
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the  .025  crystals  is  1.4/1. 7  or  .82.  This  suggest  that  the  section  from  each  boule  is  the 
same. 

And  for  Cr  we  have 

(a)  The  ratio  of  dopants  for  the  two  samples  is  .01 5/. 025  or  J> 

(b)  The  ratio  of  response  times  for  the  two  crystals  is  .39/.5 15  or  J5 

(c)  The  ratio  of  absorption,  aL,  for  the  two  crystals  is  .73/1 .2  or  Jj 

It  would  appear  that  both  NA  and  ND  change  in  proportion  to  the  Cr  concentration. 
This  is  consistent  with  our  previous  conclusion  that  the  Cr  concentration  alone  determined 
Nd  while  both  Cr  and  Mo  concentrations  determine  NA.  There  is,  however,  a  lower 
absorption  the  .015  Cr  and  .025  Cr  crystals  than  the  previous  data  for  the  double  doped 
.005  and  .01  Cr  crystals  suggested.  This  could  easily  be  due  to  the  difference  in 
absorption  through  the  boule.  It  would  be  good  to  always  look  at  crystals  taken  from 
nearly  the  same  location  in  the  boule  if  this  is  not  already  the  case. 

Conclusion 

There  is  substantial  evidence  that  in  the  double  doped  Cr-Mo  crystals  that  Cr  plays 
the  role  of  ND  the  donor,  while  both  Cr  and  Mo  play  the  role  of  NA ,  the  acceptor.  While 
the  data  is  not  100%  in  support  of  this  conclusion  the  evidence  would  warrant  a  additional 
measurements.  This  is  the  first  time,  to  our  knowledge,  that  control  on  NA  and  ND  seems 
eminent. 

The  comparison  of  our  data  on  the  new  crystals  18, 19, 20,  &21  gives  the 
following  bottom  line. 

1 .  Cr  plays  the  role  of  ND.  Only  Cr  determines  the  absorption  in  the  samples  I  have 
looked  at. 

2.  Both  Cr  and  Mo  play  the  role  of  Na.  both  Cr  and  Mo  determine  the  response  time. 

3 .  The  data  supports  the  conclusion  that  when  you  compare  the  double  dopant  to  the 
single  dopant,  the  absorption  is  given  directly  by  the  absorption  of  Cr,  i.e. 
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4 .  The  response  time  for  the  double  dopant  can  be  determined  from  the  single  dopant 

using  the  empirical  relation 

1  _  1  1 
T CrMo  TCr  ^  Mo 

This  last  result  works  extremely  well  for  the  two  cases  I  have  looked  at.  That  is, 
for 

(Cr  .005,  Mo  .005),  (Cr  .01,  Mo  .01)  and  (Cr  .005),  Mo  (.005),  and  for  (Cr  .01), 
Mo  (.01). 
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X. 


SEEDING  THE  BEAM  FAN 


In  order  to  test  the  idea  that  the  initial  scattered  light  is  weak  and  the  cause  of  a 
slow’  start  for  the  beam  fanning  phenomenon,  we  investigated  the  possibility  of  seeding 
the  fanning  with  a  more  intense  beam.  Figure  10-1  depicts  the  seeding  concept.  A 
grating  is  placed  in  front  of  the  photorefractive  crystal.  The  incident  laser  beam  now 
passes  through  the  grating  before  entering  the  crystal.  The  incident  beam  on  the  grating 
is  diffracted  into  several  beams  which  propagate  through  the  crystal  while  overlapping 
each  other.  Because  of  the  overlap,  gratings  are  written  and  energy  is  exchanged.  If 
the  grating  chosen  diffracts  more  light  than  normally  found  scattered  in  the  non-grating 
case,  the  crystal  can  be  said  to  be  seeded.  The  data  shown  in  Figure  10-2a,b  is  for  a 
series  of  gratings  of  1%  diffraction  efficiency.  The  observed  enhancement  in  speed  is  a 
modest  factor  of  3.  On  the  other  hand,  with  a  diffraction  efficiency  of  5%,  Figure  2c, 
d,  the  enhancement  in  speed  is  more  like  8  or  9.  that  is,  it  appears  that  by  providing  the 
crystal  with  a  seed  beam  of  greater  intensity  than  that  normally  found  via  scattering  by 
crystal  imperfections,  the  response  time  for  the  optical  limiting  behavior  of  the  crystal  is 
dramatically  shortened  and  brought  more  in  line  with  that  observed  for  two-beam 
coupling  experiments. 

In  order  to  investigate  the  seeding  effect  more  thoroughly  the  apparatus  of 
Figure  10-3  was  employed.  In  particular,  instead  of  using  a  series  of  gratings  to 
provide  the  seed  beam,  a  weak  seed  beam  is  intentionally  crossed  in  the  crystal  with  the 
incident  beam.  The  intensity  of  the  seed  beam  was  then  controlled  using  simple  neutral 
density  filters.  The  crossing  angle  was  controlled  by  changing  the  arrangement 
slightly.  As  seen  in  Figure  10-4a,  for  an  input  seed  intensity  of  1%  no  significant 
enhancement  was  observed.  On  the  other  hand,  for  seed  intensities  of  5,  10,  and  20% 
an  increasing  enhancement  was  noted  as  seen  in  Figure  10-4a,  b,  c.  In  addition  to  an 
enhancement  in  speed  it  is  also  interesting  to  note  that  the  transmission  was  nearly 
always  more  limited  with  the  seed  beam  than  without  (Figure  10-50. 
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DISCUSSION 


Energy  depletion  from  one  beam  into  the  other  beam  in  two-beam  coupling 

experiments  is  described  for  I01 »  I02  by  th  expression 

,  I,(0)(l  +  m) 

ll  ;“[l  +  mexp(rz)] 

r  m  I2(°Xl+m)exp(rz)  10  1 

2  [l +  mexp(rz)] 


where  1(0)  is  the  incident  intenstiy  for  beam  1  or  2,  I(z)  is  the  intensity  of  beam  1  or  2 
after  propagation  of  a  distance  z  in  the  crystal,  T  is  the  two- beam  coupling  gain 
coefficient,  and  m  is  the  ratio 

In 


m  = 


*02 


I 


01 


10-2 


The  gain  coefficient  for  zero  applied  field  is  given  by 

r=ro[l-exp(-t/x)] 

where 

27treffkBTkg  cos(20) 


ro  =- 


qXncos© 


k 

1  H - - 

k2 


10-3 


10-4 


In  this  expression  for  ro,  q  is  the  charge  of  thecarrier,  rcff  is  the  effective  electro-optic 
coefficient,  kg  is  the  grating  wave  vector,  A.  is  the  laser  light  wavelength,  0  is  the  half¬ 
angle  between  the  two  crossing  laser  beams,  and  k„  is  the  Debye  screening  wave  vector 
given  by 


K 


q2N 


eff 


kBTe 


10-5 


In  the  expression  for  k2 ,  Neff  is  the  effective  charge  density,  e  is  the  d.c.  dielectric 


constant,  and  kBT  is  the  thermal  energy. 
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The  response  time  for  he  formation  of  the  two-beam  coupling  grating  is  given 


by  the  expression 

-1  f^coYl+k;1~YYNA  |  k^T"! 

IT  ^  ctq  J_  k0  _  ^  ft  q 


10-6 


wherey  is  the  recombination  rate  coefficient,  |i  is  thhe  carrier  mobility,  NA  is  the 
acceptor  charge  density,  a  is  the  absorption  coefficient,  and  IT  is  the  total  initial 
intensity. 

As  can  be  seen  from  expression  (1)  the  change  in  intensity  is  dependent  on  the 
modulan  index  m.  Since  seeding  can  result  in  a  higher  modulation  index  than  is 
normally  present  for  beam  fanning  it  can  therefore  produce  a  greater  change  in  intensity 
via  two  beam  coupling.  This  may  explain  the  improved  or  higher  depletion  shown  in 
Figure  10-5  when  the  beam  fan  is  seeded. 

Expression  10-3  describes  the  time  behavior  of  two-beam  energy  coupling. 
The  characteristic  response  time- is  T.  As  can  be  seen  from  expression  10-6,  the 
response  time  has  an  intensity  dependence  which  is  inverse  to  the  total  intensity.  Since 
seeding  does  not  change  the  total  intensity  significantly  a  dramatic  Change  in  the  beam 
fan  time  is  not  expected.  For  example,  a  10%  or  20%  seed  would  not  be  expected  to 
alter  the  response  time  by  factors  of  3  to  8.  For  this  reason,  the  explanation  for  the 
improved  response  time  most  likely  lies  in  the  fact  that  the  formation  of  the  fan  is 
through  to  evolve  from  gratings  current.  The  addition  of  a  seed  beam  of  superior 
intensity  perhaps  encourages  the  energy  exchange  to  take  place  more  by  the  rules  of 
two-beam  coupling  than  by  the  more  complex  rules  associated  with  beam  fanning. 
Conclusions 

We  have  demonstrated  that  the  speed  of  response  for  a  beam  fanning  optical 
limiter  can  be  enhanced  by  using  a  grating  in  front  of  the  photorefractive  crystal  in  order 
to  seed  the  fanning  effect.  One  possible  extension  of  the  work  is  to  investigate  the 

possibility  of  fabricating  a  grating  right  on  the  crystal  entrance  face. 
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XI.  SUMMARY 


Our  investigations  under  this  contract  have  proved  to  be  very  productive.  We 

have: 

1 .  Developed  a  limiter,  using  a  cylindrical  lens,  with  a  103  enhancement  in 
speed.  Using  a  5  cm  x  5  cm  crystal  this  means  that  if  the  incident 
intensity  is  lW/cm2  which  normally  produces  a  0.1  second  response 
time,  the  cylindrical  technique  gives  a  100  microsecond  response  time. 
This  can  now  be  accomplished  with  an  O.D.  of  2  to  3 . 

2.  We  have  investigated  the  effect  of  a  D.C.  electronic  field  on  a 
photorefractive  crystal.  Using  an  applied  field  of  about  6KV/cm  we 
find  an  enhancement  of  a  factor  of  10  can  be  achieved.  When  coupled 
with  the  result  of  using  a  cylindrical  lens  the  combination  would  give  a 
10  microsecond  response  time  with  an  O.D.  of  2  to  3. 

3.  We  discovered  that  A.C.  applied  fields  can  cause  speed  enhancements 
on  photorefractive  SBN  and  BSKNN.  We  have,  in  fact,  developed  a 
physical  understanding  of  the  effect  of  A.C.  fields  in  these  crystals. 
The  advantage  of  A.C.  fields  is  that  the  screening  effect  accompanying 
D.C.  fields  is  eliminated.  We  have  both  computer  simulations  and 
experimental  data  to  support  our  understanding.  The  A.C.  field  can  be 
used  in  combination  with  a  spherical  focusing  lens  system.  The  two 
together  can  yield  an  enhanced  speed  of  106.  This  means  that  for  our 
incident  intensity  of  1  W/cm2,  on  the  5  cm  x  5  cm  crystal,  the  0. 1  second 
response  time  can  be  reduced  to  0. 1  microseconds. 

4.  We  developed  a  new  technique  to  measure  NA.  In  this  technique,  two- 
beam  coupling  is  used  while  we  measure  the  phase  angle  between  the 
intensity  pattern  and  the  space  change  field  as  a  function  of  applied  field. 
Our  results  were  compared  with  the  traditional  two-beam  coupling 
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technique  of  measuring  coupling  versus  crossing  angle.  The 
comparison  showed  that  both  techniques  give  the  same  answer  for  NA. 
However,  our  new  technique  is  simpler  and  more  accurate. 

5.  We  examined  writing  versus  reading  times  in  photorefractive  crystals. 
We  found  in  our  experiments  that  these  times  can  be  different.  We  also 
demonstrated  that  they  are  expected  to  be  different  using  our 
photorefractive  model.  Papers  which  state  that  these  two  times  are  the 
same  are  simply  wrong.  Similar  results  were  found  for  ordinary'  and 
extra  ordinary  polarization. 

6.  We  examined  the  response  time  of  our  crystals  versus  intensity 
Interestingly,  we  found  a  strict  1/1  dependence  on  the  response  time 
when  the  data  was  properly  taken.  This  was  true  for  several  SBN 
crystals.  BaT;03  crystals  will  soon  be  tried  as  we  showed  a  less  than  1/1 
dependence  indicating  the  presence  of  shallow  traps. 

7.  We  discovered  and  have  characterized  self  defocusing  and  focusing  by 
photorefractives.  The  defocusing  effect  in  SBN  is  very  strong,  raising 
the  divergence  of  the  incident  beam  by  3-fold  for  a  1KV  applied 
potential.  This  could  be  the  basis  of  a  “new  limiter”. 

8.  We  have  developed  a  limiting  concept  based  on  the  quadradic  electro¬ 
optic  effect.  In  this  case  our  results  demonstrate  a  last  limiter  without 
the  D.C.  screening  effect.  More  data  should  be  taken  before  these 
results  are  reported.  Here  again,  a  new  limiting  idea  was  discovered. 

9.  We  have  established  that  double  dopant  crystals  do  provide  control  over 
Na  and  Nd.  Our  evidence  is  strong  that  for  Cr,  Mo  doped  crystals,  Cr 
alone  plays  the  role  of  ND  while  both  Cr  and  Mo  play  the  role  of  NA. 
Although  the  result  is  not  the  ideal  case  of  one  dopant  acting  as  NA  and 
one  as  ND,  it  is  true  that  our  results  suggest  that  NA  can  be  varied 
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independently  of  ND  and  that  by  changing  only  the  Mo  concentration, 
Na  can  be  forced  to  change  while  holding  ND  fixed. 

1 0.  We  have  examined  the  role  of  temperature  on  the  photorefractive 
response  time.  It  is  clear  that  both  the  O.D.  and  the  reponse  time  are 
enhanced  with  higher  temperature.  More  data  must  be  taken  using  a 
combination  of  focusing,  field,  and  temperature  before  a  report  is  made. 

1 1.  We  have  developed  a  seeding  method  to  increase  the  photorefractive 
response  time  by  a  factor  of  5  tolO  without  loss  of  O.D. 

We  believe  that  the  evidence  is  strong  that  an  effective  photorefractive 
limiter  is  NOW.  The  next  step  is  to  try  a  device  for  each  “concept”  investigated  or 
demonstrated,  and  to  compare  with  the  needed  response  time  and  O.D.  I  would 

be  happy  to  Support  any  effort  by  the  ARMY  to  construct  a  photorefractive  optical 
limiter  based  on  the  results  of  our  work. 

Finally,  we  have  three  strong  suggestions  for  device  fabrication  based  on  our 
studies.  The  first  is  to  use  a  cylindrical  lens,  spherical  lens  combination.  Consider 
an  input  beam  5  cm  x  5  cm.  A  spherical  lens  system  reduces  the  beam  to  .5  cm  x 
.5  cm,  while  a  cylindrical  lens  system  takes  the  beam  to  .5  on  x  50  microns.  The 
total  reduction  in  area  is  104.  As  a  result  the  response  time  is  104  times  faster!  If 
we  now  apply  a  D.C.  electric  field,  we  can  produce  a  105  enhancement  in  speed. 
This  device  would  not  suffer  from  screening  and  would  have  a  2  to  3  O.D.  This 
works  TODAY,  limiting  at  6  uj/cm2  with  an  O.D.  of  3  to  4.  In  addition,  seeding 
the  beam  fan  could,  in  principle,  bring  improvement  higher. 

The  second  device  is  based  on  using  a  spherical  lens  system  coupled  with  an 
A.C.  electric  field.  Using  the  same  parameters  outlined  above  and  a  BSKNN 
crystal  and  10*  enhancement  in  speed  is  possible  with  an  O.D.  of  2  to  3.  While  this 
system  is  simpler  than  the  cylindrical  lens,  the  laboratory  results  were  never  as 
clean  cut.  However,  it  is  promising  and  should  be  tried. 
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The  third  device  is  based  on  the  quadradic  electro-optics  effect.  We  found 
exciting  results  using  an  A.C.  square  wave  on  a  SBN  sample  with  La  doping.  At  room 
temperature,  this  sample  was  above  but  near  the  Curie  temperature.  The  experimental 
results  showed  a  modulated  scattering  or  fanning  phenomena  that  followed  the  temporal 
behavior  of  the  driving  field.  This  system  protected  effectively  with  spherical  lens, 
except  during  the  zero  applied  field  periods  of  the  cycle.  However,  a  combination  of 
two  crystals,  each  at  a  different  driving  frequency,  would  illuminate  this  problem. 
Here,  again  we  found  a  limiter  with  promise  for  today. 

In  addition  to  the  H  major  results  of  our  work  it  should  be  noted  that  at  least  22 
publications  in  major  journals,  20  publications  in  proceedings,  and  about  40  talks  at 
conferences  were  also  among  the  achievements  of  the  grant.  Moreover,  additional 
papers  will  be  submitted  over  the  next  six  months. 

It  should  be  noted  that  nearly  all  my  work  was  carried  out  with  the  help  of  many 
of  the  scientists  at  the  Night  Vision  Electro-optics  Center  part  of  which  is  now  the 
Army  Research  Laboratory.  I  am  grateful  for  the  support,  collaboration,  and  friendship 
of  Ed  Sharp,  Gary  Wood,  Bill  Clark,  Mary  Miller,  and  Byong  Ahn.  Of  course  none 
of  my  work  would  have  been  possible  without  the  crystals,  scientific  expertise,  and 
collaboration  of  R2  Neurgaonkar.  Also,  I  would  like  to  add  a  special  note  of  thanks  to 
Wolfgang  Elser  and  L.  N.  Durasula,  who  supported  this  work  from  beginning  to  end 
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Appendix  I :  Derivation  of  the  Space-Charge  Field 


In  this  section,  the  steady-state  and  time-dependent  solutions  for  the  space-charge  field 
will  be  derived  from  the  linearized  set  of  equations  given  by  (2.6.a)-(2.6.h): 


dK 


Jt  =sI0(Nd-N+0)-  yRne0N+0 


(I.l.a) 


dNt 


Jo  =  QeVene0E0  - 


h  =  (le\ie(<ne0E\  +  ne\E0)  +  i\iekBTkgne\  > 


(I.l.b) 

(1.1. C) 

(I.l.d) 


0 


d_K  _dne o 

dt  dt 


(I-l.e) 


..  •  __  (dN+i  dnel) 

V 1  q\dt  dt 


0  =  N0~ne0  ~NA 


ik8E  i  =  i) 


d.i.g) 

(Ilh) 


First,  an  expression  for  nei  must  be  obtained  which  eliminates  Combining  the  time 
derivative  of  (I.l.h), 


dE  i 

8dt  e  \dt  dt  j 


with  (I.l.d)  and  (I.l.f),  one  obtains 


(1.2) 
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(1-3) 


=  qAi 

dt  8 


( ne0El  +  nelE0 ) 


i[lekBTkgnel 


Solving  this  for  ne  j  yields 


ne\  = 


ik  edE, 

— - — y  +  m  k.n.nE, 
qe  dt  g  e 0  1 

J _ L 

TZ>  T£ 


(1.4) 


where  the  time  constants  in  (1.4)  are  defined  as 


xn  =  mean  diffusion  time  =  - -  , 

L*  . .  i  m  i  2 


(I.5.a) 


x£  =  mean  drift  time  =  — t-tt  .  (1.5  .b) 

Assuming  neg  to  be  constant,  the  time  derivative  of  (1.4)  is  given  by 


ikged2El 
Qe  dt 2 


dE{ 

+  ''MAo* 


(1.6) 


At  this  point,  (I.l.b)  and  (1.6)  can  be  substituted  into  (1.2)  to  obtain  an  expression  involv¬ 
ing  Afj+and  ne\  and  time  derivatives  of  E\: 


smtQ (,Nd-M*0)  -  sl0 - yR  (ne0V\  +ne[N*0) 


ik.ed-E . 

— - -*  +  /H  *„ n 

9,  d,2  '  g 


dE{ 

e0Jt 


ikgzdEl 
Qe  dt 


(1.7) 


Using  (I.l.h)  and  (1.4)  to  eliminate  Nj+  and  nei  from  (1.7)  results  in  the  following  second- 
order  differential  equation  for  E\. 
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d~E,  111  ;  dE}  1  1  ;  Q  1  ; 

tt  +(r+r+rT)7r  =  /n(rT)m/0(^*-^0“^) 

«*  +  di  lD  l£  V  Vi  I  D  X/X£  ,A^E  X£>  X£ 


,  d.8) 


where  the  time  constants  in  (1.8)  are 


%di  =  dielectric  relaxation  time  = 


qeVeneO 


(I.8.a) 


%j  =  excitation  time  = 


*4) +  ^Rne  0 


(I-8.b) 


x.  = 


+  sI0  +  2yRne0  +  y rNa 


(I.8.C) 


Since  neQ  is  constant,  (I. l.a)  and  (I.l.e)  can  be  used  to  obtain  the  following: 


dNt  dn„ 


sI0(ND-N*0)-yRne0Nl=Ti°  =^'°  =  0. 


(1.9) 


Assuming  NA  »  ng0  for  low  irradiances,  the  resulting  expression  for  neg  is 


ne0  ~ 


sIo  (nd~na) 

yR»A 


(1.10) 


At  equilibrium,  the  time  derivatives  in  (1.8)  vanish.  Hence,  the  steady-state  expression  for 
E\  is  given  by 


E\s 


Jj^(^--^-)s»il0(ND-ne0-NA) 
&  D  E 


(^U+  1 


(1.11) 


\Xdi  xixd  xixe 


) 


For  low  irradiances  ( NA  »  ng0  and  yRNA  »  sIQ),  (1.1 1)  can  be  simplified  to 


E\s 


=  mE 


N 


ED+iE0 


Eq  +  i  (Ed  +  EN )  J 
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(1.12) 


where  Ed  and  are  defined  as 


,  _tBn, 
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4e 


r  Qe^eff  Qe^A 

£*  =  ^r  =  ij 


na  v> 


y+ND-NAJ 


Note  that  (1.12)  can  be  written  in  the  form 


Ef  =  mEj ♦  , 

where  the  amplitude,  Esc,  and  the  phase  shift,  <{),  are  given  by 


„  «.  ,  4+4 

""  1  (£«+£d)2  +  4’ 


tan(|)  = 


Ed f 
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En  EdEN/ 


(LI  3. a) 


(I.13.b) 


(1.14) 


(I.15.a) 


(I.15.b) 


The  transient  character  of  E\  is  found  by  solving  the  homogeneous  form  of  (1.8).  This  can 
be  written  as 


d2E{  i  dEx 
1?  + 


where  Tj  and  I 2  are  defined  by 


(1.16) 


(I.17.a) 


T  T  .. 
+  di 


XIXD 


XIXE 


(I.17.b) 
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The  complex  quantities  xa  and  xb  can  be  defined  such  that 

1  1  1 
=  —  +  — 


1  1  1 


The  solution  to  (1.16)  then  takes  the  form 


(I.18.a) 


(1.1 8.b) 


£j(r)  =  E\s  +  Ae  t/X°  +  Be  ,A\ 


(1.19) 


where  A  and  B  are  constants  to  be  determined  by  initial  conditions.  By  assuming  x+  «  %di, 
x+  «  xD  and  x+  «  x/,  (I.18.a)-(I.18.b)  can  be  solved  for  xa  and  xb  in  terms  of  T\  and  T2: 


(I.20.a) 
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x 


b 


(I.20.b) 


Since  Re{ l/xfl)  «  Re  (l/xb) .  the  term  in  (1.19)  involving  xb  can  then  be  neglected 
since  it  decays  much  faster  than  the  term  involving  xa.  Thus,  (1.19)  becomes 


£j  (0  =  E\s  +  Ae 


(1.21) 


The  real  quantities  x  and  to  can  now  be  defined  such  that 


1 


x 


a 


1  . 

-  +  zco  . 
X 


Using  (I.17.a)-(I.17.b),  the  following  expressions  for  x  and  to  are  obtained: 


(1.22) 
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This  initial  value  problem  is  completed  by  determining  the  constant  A  in  (1.21).  If  no  grat¬ 
ings  exist  in  the  photorefractive  material  at  time  t  =  0,  then  the  initial  condition  is 
£j  (0)  =  0,  so  the  time-dependent  form  of  the  space-charge  field  is  given  by 


Ex(t)  =  E\s  [  1  -  e~t/xe~im]  , 


(1.24) 


where  x  and  to  are  given  by  (I.23.a)-(I.23.b),  respectively. 
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APPENDIX  I-B 


DERIVATION  OF  EQUATION  ( 4-8  ) 

The  two  waves  R  and  S  (reference  and  signal  waves  respectively  )  intersect  in  the 
volume  of  the  real  time  holographic  medium.  The  amplitudes  of  two  waves  are 
deduced  as 

dR/dz  =  -XoS  (1-1) 

dS/dz  =  XoR  d-2) 

where  x0=  rmtAnsin(p/(X.cos0)  with  m  =  2RS/(R2  +  S2)  is  the  corresponding  beam 
coupling  constant.  R(0)  =  R0,  S(0)  =  S0  are  the  boundary  conditions. 

By  introducing  r=47tAnsin(p/(Xcos9)  and  substituting  the  expression  of  m  into  eq.  ( 
I- 1 )  and  ( 1-2  )  we  obtain 

dR/dz  =  -(l/2)rRS2/(R2  +  S2)  ( 1-3  ) 

dS/dz  =  (l/2)rR2S/(R2  +  S2)  ( 1-4 ) 

We  have  IR  =  R2,  Is  =  S2  and  R2  +  S2  =  R02  +  S02. 

Therefore,  the  R  and  S  amplitudes  satisfy  to  the  equations: 

dR/dz  =  -( 1/2)TR  +  (D2)R3/(R02  +  S02)  ( 1-5  ) 

dS/dz  =  (l/2)rs  -  (r/2)S3/(R02  +  S02)  ( 1-6 ) 

The  waves  intensities  emerging  from  the  material  ( z=Lbff )  and  solutions  of  eq.  (1-5) 
and  (1-6)  are  given  by 

Is  =  Iso(  V  +  S02 )/[  R02  +  S02exp(rLeff)  ]  ( 1-7  ) 
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Ir  =  Iro(  Ro2  +  So2 )/[  Ro2  +  S02exp(TLeff)  ]  ( 1-8  ) 

If  Ilc  and  i!  are  transmitted  intensities  of  the  signal  beam  with  and  without  coupling, 
respectively,  then  signal  beam  (  suppose  it  is  a  weak  beam  ),  Ij  will  experience  an 
increase  in  intensity  given  by 

I]c/Ij  =  ( Is  with  reference  beam )/( Is  without  reference  beam ) 

=  ( Ro2  So2  )exp(FLeff)/[  Rq2  +  SQ2exp(rLef-f)  ] 

=  ( 1  +  c  )exp(TLeff)/[  1  +  cexp(TLefr)  ] 
with  c  =  S02/R02.  This  is  eq.  ( 4-8  )  used  in  Chapter  4. 
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APPENDIX  II 

ANALYTICAL  SOLUTION  FOR  SPACE-CHARGE  FIELD  WITH 
AN  APPLIED  SQUARE  WAVE  AC  FIELD 


Using 


d_Ei  1  dE,  1 

at2  +  T,at  +  t2’e'-c 


1/T,=A1+  iAjEq 
1  /T22  =  B,  -*-iB2E0 
c  =  i(C,+iC2E0)  =  -C2E0  +  iC, 

where  A,,  A2,  B,,  B2,  CJf  and  Cj  are  real  numbers. 


Then 


-V(T22)  =  -(B,+iB2E0)/(A,+  iA^) 

=  -  ( ( A,B,  +  A2B25d2)  +  iE0(  -A2B,  +  A,B2 )  ]  / ( A,2  +  A22^,2 ) 
=  -1/t  -  ico 

T22c  =  (-C2E0  +  iC1)/(B,+iB2E0) 

=  ( ( *  BjCj  +  C,B2)  Eo  +  i  ( B,C,  +  B2Cj^2)  ] /( B,2  +  B22Eo2 ) 

=  ESc(  cos  9  +  i  sin  <p ) 


=  Etce» 


9  -  tan-1  {  ( B,C,  4  B2C2E02 )  /  [  ( -  BlC2  4  C,B2  )E0]) 


(ii-D 


( n-2 ) 
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Introduce  a  non-dimension  physical  quantity  E,  ( t )  as 

E,(t)  =  EscE,'(t) 

If  we  apply  a  square-wave  AC  field  across  the  crystal  the  amplitude  and  period  of 

which  are  E0  and  T,  respectively,  the  E,  in  the  end  of  the  first  half  period  would  be 
E,  ( t  =  T/2 )  =  - e *e -T/(2t>e '•»T/2+e  (  n-3 ) 

During  2nd  half  period  ( T/2  <  t  £  T ),  the  applied  field  is  -  So  now 

T2’c  =  [  ( -B,Cj  +  C,B2 )  ( -  E„ )  +  i  ( B,C,  +  BjCjE05  )  /  ( B,2  +  B2V> 

=  Ek(  -cosep  +  i  sirup ) 

=  -Ekc-"> 

*  T,  /  ( Tj  )  =  -  (  Bj  -  i  B2Eo  )  /  (  A,  -  i  A2E0 ) 

=  -  [  ( A,B,  +  A2B2Eo2  )  -  >  H0(  -A2B,  +  A,B2 )  ]  /  ( A,2  +  a/e/  ) 

=  -1  /t  +  ito 
Then  the  solution  is 

E/( t  )  =  I E,‘ ( t  =  T/2 )  +  e‘,v]  •  e ^'-T/2)/T  •  e  to(,  T/2>- e1’ 

=  l -ekpeT/(2t,e,WT'2  +  e,<?  +  e* ]  •  e*(,-T/2)/T .  eto<« 'T/2> . e*  ( T/2  <  t  <  T ) 

Thus 

E,‘  ( t  =T )  =  -  c‘ VT/t  +  e*  ■  eT/  ( *  >  •  ctoT'2  +  e* .  e-T/2t .  cu»t/2  ^ 

By  using  the  same  procedure  as  above,  we  can  get 

E,(  t  =  3T/2 )  =  /e,T/(2,)ctoT/J  +  e*  eTn  +  c‘v  e'T/T  -  e-i?e-T/(2,)ctoT/2 

-  e*^e'T/2t  ei(*V2  +  e1’ 
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E/(t  =  2T)  =  -ei<p e‘2T/ x  +  eV37'* eiwT/2  +  eicp e37/2x eicoT/2  - ei<pe*T/x - ei<peT/T 


+  ei<p  e’T/  2x  e‘“T/2  +  e  i(pe  T/2tetoT/2  -  e‘ilp 


E,’  ( t  =  5T/2  )  =  -ei<p  e'57/2xe'ia)7/2  +  e’V17*  +  e'i(p  -  e^e-3T/2xeicoT/2 

.  ei(Pe-3T/2te'i“T/2  +  eilPe-2T/2T  +  e-^e-2T/2r  _  g-i<p g-T/ 2t  e-ia)T/2 


ei(pe-T/2te.i0)T/2  +  ei<P 


Now  we  want  to  find  the  recurrence  in  the  expression  of  E[  (  nT  )  (  n  is  a  positive 
integer ).  Look  at  E,’  ( 2T  )  and  rewrite  the  expression  as 


Ej’  (t  =  2T)  =  -e^e™ ( 1  +  e‘7/x )  +  ei(p eiCl)T/2 e'7/2x ( 1  +  e'7/x )  - ei<p ( 1  +  e‘7/x ) 


+  e‘i<p  •  eimT/2  •  e‘7/2x  ( 1  +  e"T/x ) 

From  this,  we  can  derive  an  expression  for  E,  ’  as  t  =  nT,  which  is 

E;  ( t  =  nT )  =  [  1  +  e'T/T  +  e 27/x  + . . +  e'(n'1)7/x]  •  ( -e1  VT/X  +  ei<f,-eiCl)T/2 

.  e'T/2x  -  e"i<p  +  e'i<p-  e'“T/2  •  e‘7/2x) 


=  (  1  -  e*nT/  x )  /  ( 1  -  e* w  T )  •  ( -eKp  e*w  x  +  e,,,,  •  e 
+  e'i<f*-e'“T/2-e'T/2x ) 


-T /T  ^  „*T/ X  ,,*9 .  JwT/2  .  e-T/2x  _  e-i9 


(n  =  0,1,2, . )  ( II-4 ) 

Next  let’s  look  at  E[  ( t  =  5/2T )  and  rewrite  it  as 

E,'  ( t  =  5T/2 )  =  e* ( 1  +  e’T/x  +  e’27/x )  - ei<p ei£oT/2  eT/2x  ( 1  +  e'T/x+  e2T/x )  + 
e-‘VT/x  •  (  1  +  e'T/t )  -  e‘i(t*  e iuT/2e‘T/2x ( 1  +  eT/x ) 

=  ( 1  +  e"T/x+  e'2T/x )  ( e*  -  ei<p  e4*7 11  e7/2x )  +  ( 1  +  e7/x ) 

( e'i<pe'T/x  -  e‘i<t>  e’’“7/2  e'7/2x ) 
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Also,  from  this,  we  can  predict  the  expression  for  E,  as  t  =  (  n  +  1/2  )  T, 
which  is 

E,'  ( t  =  ( n  +  1/2 )  T  )  =  [  1  +  e‘T/x  +  e2T/x  + . +  e'nT/x]  -(ei<p- 

ei<Pe,C0T/2.e-T/2x)  +  [1+e.T/x  +  e-2T/t+ . + 

e-(n-l)T/x]<e-iVT/x_e-i<pe-ia,T/2e-T/2x) 

=  ( 1  -  e'(n+1)T/x )  / ( 1  -  e'T/T )  •  ( ei<f -  ei<p e'i£l)T/2  eT/2x ) 

+  ( 1  -  e'nT/T )  /  ( 1  -  e'T/t )  •  ( e'iipe'T/T  -  e'i(p  eifflT/2  e‘T/2x ) 

(n  =  0, 1,2, . )  ( H-5  ) 

Eq.  ( II-4  )  &  ( II-5  )  are  the  basic  equations  (5-1  )  &  (  5-2  )  in  Chapter  5. 
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APPENDIX  IH 


DERIVATION  OF  EQUATION  FOR  PHASE  SHIFT  MEASUREMENT 

Two  laser  beams  with  ordinary  polarization  and  equal  intensity  I0  intersect  inside  a 
crystal,  which  results  in  an  energy  exchange  signal  I£and  a  diffraction  signal  Id. 

The  diffraction  efficiency  r|  is  given  by  ( Ref.  [1],  [2] ) 

n  =  Id/I0  =  sin2(Kd)  ( III- 1 ) 

where  k  is  the  coupling  constant  and  d  is  the  thickness  of  the  crystal. 

According  to  Ref. [7],  both  the  real  and  imaginary  part  of  the  space-charge  field 
contribute  to  the  diffraction  efficiency.  Therefore,  Kin  eq.  (III-l)  is  given  by  (  Ref. 
[1],  [2],  [7] ) 

K  =  7t-An/Xcos0  - 

=  (7t/2?i  cos0)(r^3  refr)[(ReEsc)2+(ImEsc)2]1/2  ( III-2  ) 

where  ReEscand  ImEscare  real  and  imaginary  part  of  the  space-charge  field, 
respectively.  If  we  write 

c  =  {%12X  cos0)(no3  refr)  ( III-3  ) 

Then  from  eq.  (III-l)  and  (III-2)  we  can  obtain 

[(ReEsc)2+(ImEsc)2]  m=  (l/c-d)[sin->(Tl m)}  ( HI-4 ) 

In  other  hand,  the  energy-exchange  efficiency  is  ( Ref.  [1] ) 

e  =  I£/I0 

The  gain  of  the  detected  beam  is  governed  by  the  equation  ( Ref.  [1]  and  Appendix  I ) 
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( Io+g/ 10  =  l+e  =  (l+l)erd/(l+erd)  ( III-5  ) 

where  T  is  the  gain  coefficient,  which  is  only  due  to  the  imaginary  part  of  the  space- 
charge  field  (  Ref.  [7] ).  T  is  given  by  ( Ref.  [4] ) 

T  =  (2nfk  cos0)(no3  reff)(ImE.c) 

By  using  eq.  (III-3)  T  becomes 

r  =  4c(ImEsc)  ( III-6  ) 

By  substituting  eq.  (III-6)  into  eq.  (III-5),  we  can  obtain 

ImEsc=  (l/4)(l/c-d)-ln[(l+e)/(l-8)]  ( III-7  ) 

The  phase  shift  cp  between  the  grating  and  the  fringe  pattern  is  determined  by  the 
equation 

tancp  =  ImEsc/ReEsc 

Therefore,  by  using  eq.  (III-4)  and  (III-7)  we  will  get 
sintp  =  tantp/(l+tan2(p)1/2 

=  ImEsc/[(ReEsc)2+(ImESc)2]I/2 

=  ln[(l+e)/(l-£)]/[4sin*I(tl1/2)]  ( III-8  ) 

For  small  coupling  which  is  the  case  in  our  experiment,  both  tj  and  e  are  much  less 
than  1.  So  mathematically,  by  using  Taylor  expansion,  we  can  get 
ln[(l+e)/(l-e)]  =  2e  fore<l 

and  sin->(ri,/2)  =  T|I/2  forricl  (III-9) 

By  substituting  eq.  (III-9)  into  eq.  (III-8 ),  we  obtain 


179 


sincp  =  e/2r|1/2 


=  (V  i0)/2(V  I0)1/2 

=  I€/2(IdI0)1/2  . 

Equation  (III- 10)  is  just  equation  (6-1)  presented  in  Chapter  6. 


(III-10) 
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APPENDIX  IV 


NUMERICAL  METHOD  AND  PHYSICAL  PARAMETERS 

Numerical  Calculation  for  the  Space-Charge  Field 

The  equation  of  space-charge  field  is  time  dependent  second-order  differential 
equation,  the  form  of  which  is 

d2Ei(  t  )/dt2  +  A  ( t )  dEi(  t )/  dt  +  B(  t )  Ei(  t )  =  C(  t ) 
where  Ej(  t ),  A(  t ),  B(  t )  and  C(  t )  are  all  complex. 

Firstly,  we  have  to  separate  the  real  and  imaginary  part  of  the  above  equation. 
Secondly,  we  must  separate  the  second-order  differential  equation  into  two  first-order 
differential  equations  by  a  change  in  variables. 

Thirdly,  we  utilize  Euler’s  method  as  a  numerical  method  for  solving  differential 
equation. 

Fourthly,  we  write  a  C  program  to  generate  the  numerical  solution. 

Physical  Parameters  Used  for  Computer  Programs 

Dielectric  constant  €  =  3000 
Temperature  T  =  300  k 
Optical  wavelength  1  =  4580  A° 
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Photoionization  cross  section  s  =  5-10*18  cm^ 
Recombination  coefficient  yR  =  1-10*9  cm^/s 

Mobility  |i  =  0.5  cm^/vs 


Full  crossing  angle  20  =  10° 


APPENDIX  V 


COMPUTER  PROGRAM  FOR  CALCULATING 
PHOTOREFRACTIVE  SPACE-CHARGE  FIELD 

/*  Program:  Space-charge  field  with  applied  square-wave  AC  field. c  */ 

#include  <stdio.h> 

#include  <math.h> 

#include  <stdlib.h> 
void  main(void) 

{ 

int  i,  points=250000,  j,  1,  m; 

long  double  constO=3.861935e+l,  Tplus=8.0e-9,  s=5.0e-19,  gammaR=2.5e-8, 

gammaRNA=1.25e+8,  mu=5.0e-l,  k=5.0e+3,  NDminusNA=5.0e+16, 
const l=1.66e+9, 10,  omega,  Tc,  TD,  hc=1.988e-15,  nO,  const2,  Tdi,  TI, 
TE,  Al,  A2,  Bl,  B2,  Cl,  C2,  lemda=4.58e+3,  I=2.4e+1,  E0=3.0e+3, 
delt=7.0e-9,  RE=0.0e+0,  IE=0.0e+0,  dRE,  dIE; 
static  long  double  ReEl[5000],  ImEl[5000]; 

I0=I*lemda/hc; 

n0=s*I0*(NDminusNA)/(gammaRNA); 

TD=constO/(mu*k*k); 

Tdi=constl/(mu*nO); 

TI=l/(s*IO+(gammaR)*nO); 

TE=  1  /  (k*mu  *E0) ; 

Tc=Tdi*((l+Tplus/TD)*(l+Tplus/TD)+(Tplus/TE)*(Tplus/TE))/ 
((l+Tplus:):Tdi/(TD*TI))*(l-i-Tplus/TD)+Tplus*Tplus*Tdi/(TE*TE*TI)); 
omega=(Tplus/TE)*(Tdi/TI- 1 )/ 

(Tdi*((l+Tplus/TD)*(l+Tplus/TD)+(Tplus/TE)*(Tplus/TE))); 

const2=s*I0*((NDminusNA)-n0)/(constl*k); 

Al=l/TD+1/Tdi+1/Tplus; 

A2=1/TE; 

Bl=l/(TD*TI)+l/(Tplus*Tdi); 

B2=1/(TI*TE); 

Cl=-const2/TE; 

C2=const2/TD; 

dIE=(C2*Al-Cl*A2)/(Al*Al+A2*A2); 
dRE=(C  1  *  A 1  +C2*A2)/(A  1  ^  A 1 + A2*  A2); 

ReEl[0]=0.0e+0; 

ImEl[0]=0.0e+0; 
for  (1=0;  1<=360;  1++) 
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{ 

for  (i=(2*l*points);  i<((2*l+l)*points);  i++) 

{ 

RE=RE+delt*(dRE); 

dRE=dRE+delt*(Cl-Al*dRE+A2*dIE-Bl*RE+B2*IE); 

IE=IE+delt*(dIE); 

dIE=dIE+delt*(C2-Al  *dIE-A2*dRE-B2*RE-B  PIE); 

} 

ReEl[2*l+l]=RE; 

ImEl[2*l+l]=IE; 

printf("\n%5d  %Le  (2*1+1),  ReEl  [2*1+1],  ImEl  [2*1+1]); 

for  (j=((2*l+l)*points);j<(2*(l+l)*points);j++) 

{ 

RE=RE+delt*(dRE); 

dRE=dRE+delt*(-C  1  -A  1  *dRE-A2*dIE-B  1  *RE-B2*IE); 
IE=IE+delt*(dIE); 

dIE=dIE+delt*(C2-Al*dIE+A2*dRE+B2*RE-B  1  *IE); 

} 

ReEl[2*l+2]=RE; 

ImEl[2*l+2]=IE; 

printf("\n%5d  %Le  %Le",  (2*1+2),  ReEl  [2*1+2],  ImEl  [2*1+2]); 
printf("\n\n"); 

} 
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